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ABSTRACT 

\ 

Theoretical  and  experimental  structural  response  investigations 
of  space -vehicle -type  structures  under  suddenly  applied  external  sur¬ 
face  loads  are  described.  The  simulation  of  a  simultaneous  impulsive 

load  by  a  traveling  load  such  as  produced  by  an  explosive  is  analyzed  j 

1 

for  the  string  and  membrane.  Three  dynamic  buckling  problems  are  ! 

investigated:  (1)  dynamic  plastic -flow  buckling  of  fiat  plates  due  to 

in-plane  flow,  (2)  dynamic  elastic  buckling  of  a  thin  cylindrical  shell 

under  axial  impact,  and  (3)  dynamic  buckling  of  cylindrical  shells  of 

a  strain-rate  sensitive  material.  A  scheme  for  correlating  the  results 

of  structural  response  investigations  concerned  with  dynamic  failure 

loads  of  structures  is  presented  and  a  brief  review  of  available  results  * 

is  given. 
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1.  INTRODUCTION 


The  work  covered  in  this  report  is  an  extension  of  that  described 
in  [  1.  l]  and  [1.2],  The  specific  motivation  for  the  investigations  is 
discussed  in  Volume  I  of  [  1 .  1  ] . 

The  work  reported  here  falls  in  three  major  parts:  simulation 
of  simultaneous  loads  by  traveling  loads,  investigations  of  dynamic 
buckling  of  structures,  and  correlation  of  structural  response  infor¬ 
mation. 

%  , 

Simulation  of  simultaneous  loads  by  traveling  loads  is  of 
considerable  interest  because  of  the  relative  ease  of  using  explosives 
to  obtain  impulsive  loads.  In  Section  2,  the  response  of  the  string 
and  the  membrane  to  traveling  loads,  as  a  function  of  load  velocity, 
is  investigated.  It  is  found  that  detonation  velocities  are  sufficiently 
high  to  give  a  good  approximation  to  a  simultaneous  load. 

Buckling  has  been  found  to  be  a  significant  response  mode  for 
space  vehicles  under  suddenly  applied  external  surface  loads.  In 
the  past  year  three  dynamic  buckling  problems  have  been  investigated. 
These  are  reported  in  Sections  3,  4,  and  5.  Section  3  describes 
dynamic  plastic -flow  buckling  of  a  plate  due  to  in-plane  flow.  Section  4 
describes  dynamic  elastic  buckling  of  a  thin  cylindrical  shell  under 
axial  impact.  Section  5  describes  dynamic  buckling  cf  a  cylindrical 
shell  of  a  strain-rate  sensitive  material. 

Correlation  of  structural  response  information  related  to 
failure  loads  of  space -vehicle -type  structures  is  desirable  to  make 


Numbers  in  brackets  indicate  the  references  listed  at  the  end  of 
each  section. 
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the  results  more  accessible  and  to  indicate  areas  where  additional 
work  is  required.  In  Section  6  it  is  shown  that  for  dynamic  loads, 
the  amplitude -impulse  plane  is  a  convenient  and  meaningful  way  to 
present  failure  loads.  Using  this  format,  a  brief  summary  of 
available  results  is  given. 

In  Section  7  some  comments  are  made  on  the  separation  of 
the  response  of  space  vehicles  to  suddenly  applied  surface  loads 
into  shock  response  and  structural  response,  and  on  prediction  of 
accelerations  transmitted  to  internal  components. 


REFERENCES 


1  Abrahamson,  G.  R.  and  A.  L.  Florence,  "Investigation  o£ 
Response  of  Simplified  ICBM-Type  Structures  to  Impulsive 
Loading  (U),"  Volumes  I  and  II,  AFSWC-TDR-62-94,  Contract 
AF  29(601)-4329,  November  1962  (SRD)  and  December  1962 
(Unclassified). 

2  Abrahamson,  G.  R. ,  A.  L.  Florence  and  H.  E.  Lindberg, 
"Investigation  of  Response  of  Simplified  ICBM-Type  Structures 
to  Impulsive  Loading  (U),"  Stanford  Research  Institute  Final 
Report,  AFWL-TDR -64-22,  Contract  AF  29(601)-4329,  July  1964. 
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2.  SIMULATION  OF  SIMULTANEOUS  IMPULSIVE  LOADS 
ON  STRINGS  AND  MEMBRANES  WITH  EXPLOSIVES 


2.  1  INTRODUCTION 

Explosives  are  a  convenient  means  of  simulating  simultaneous 
impulsive  loads  on  structures.  The  explosive  imparts  an  impulse  by 
means  of  a  high-pressure  pulse  acting  within  a  narrow  reaction  zone 
at  the  detonation  front  which  travels  over  the  surface  away  from  the 
initiation  point  at  the  detonation  velocity.  Consequently,  elements 
of  the  structure  receive  impulsive  velocities  successively  and  not 
simultaneously  as  desired.  Good  simulation  of  a  simultaneously 
applied  impulse  over  an  area  can  be  expected  if  the  detonation  velocity 
is  sufficiently  large,  but  some  criterion  is  required  to  decide  when 
this  velocity  is  large  enough. 


In  this  section  the  responses  of  a  stretched  infinite  string  and 
a  stretched  infinite  membrane  to  traveling  concentrated  forces  repre¬ 
senting  a  detonation  front  are  investigated.  The  string  has  a  concen¬ 
trated  force  of  magnitude  2P  suddenly  applied  to  it  which  is  supposed 
to  separate  immediately  into  two  equal  forces  of  magnitude  P  and 
travel  in  opposite  directions  each  at  velocity  V  (Fig.  2, 1).  The 

membrane  has  applied  to  it  a 


a  a 


FIG.  2.1  STRING  PROBLEM 


circular  ring  load  of  magnitude 
P  per  unit  length  the  radius  of 
which  increases  at  a  constant 
velocity  V  (Fig.  2.2).  Of 
particular  interest  are  the  dis¬ 
placement  and  velocity  distri¬ 
butions  imparted  to  the  string  and 


membrane  when  the  loads  are  moving  supersonically  relative  to  the 


wave  velocities.  The  velocity  distributions  are  compared  with  the 


constant  velocity  distribution 
resulting  from  the  whole  impulse 
being  applied  uniformly  over  the 
length  or  area  traversed  by  the 
load. 

Past  studies  originated  with 

a  treatment  by  Goodier  [2.1]  of 

a  stretched  semi-infinite  string 

subjected  to  a  concentrated  force 
FIG.  2.2  MEMBRANE  PROBLEM  ,  .  , 

which  runs  on  to  it  at  the  support 

and  travels  at  a  constant  velocity. 

It  was  shown  that  the  velocity  distribution  approaches  that  due  to  an 
ideal  impulse  covering  the  portion  of  the  string  traversed  by  the  load 
as  the  velocity  of  the  force  becomes  large.  Florence  [2.2]  treated 
the  corresponding  problem  for  the  beam  and  included  the  problem 
corresponding  to  the  string  problem  studied  here.  A  similar  con¬ 
clusion  was  drawn  but  additionally  the  velocity  distribution  was  found 
much  more  uniform  when  initiation  was  away  from  a  fixed  support. 
This  permitted  very  uniform  distributions  of  velocity  with  practical 
values  of  beam  wave  velocities  and  detonation  velocities.  Also,  it 
was  found  desirable  to  have  a  detonation  velocity  supersonic  relative 
to  the  beam  wave  velocities. 

It  is  shown  here  that  with  central  initiation  and  with  practical 
values  of  wave  and  detonation  velocities,  a  uniform  velocity  distri¬ 
bution  on  the  string  and  membrane  is  achieved.  Also,  results  are 
presented  to  describe  the  distributions. 

2.2  STRING  THEORY 

Let  the  two  constant  loads  P  originate  at  time  t  *  0  at  the 

origin  0(x,  y)  (Fig.  2. 1)  and  separate,  each  with  a  constant  velocity 
V  .  Let  the  string  have  mass  m  per  unit  length,  be  stretched  by  a 
force  s  and  be  at  rest  at  time  t^O.  Then  the  differential  equation 


of  motion,  the  initial  conditions,  and  boundary  conditions  are 


6(x  -  Vt)  x  >  0  (2. 1) 

dx  c  St  me 

y(x,0)  =  |f  (x,0)  =  0  (2.2) 

and 

|£(0,t)  =  0  Lim  y(x,  t)  =  0  (2.3) 

x  -*® 

1/2 

where  c  =  (s/m)  is  the  string  wave  velocity,  and  6  is  the  Dirac 
delta  function. 

Applying  the  Laplace  transformation  to  (2. 1)  and  (2.  3),  and 
making  use  of  (2.2)  gives 

-  (p/c)2y 
dx 

(0.  p)  =  0 


P  -px/V 
—  •  —  ■  ■  c 

me  V 


(2.4) 


Lim  y(x,  p)  =  0 

x  ■*  • 


(2.5) 


The  solution  of  (2.  4)  satisfying  (2.  5)  is 


y(*.  P) 


P 

m 


1 


P 


(2.6) 
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Inverting  the  transform  (2.  6),  noting  that  the  only  singularity 
is  a  simple  pole  at  the  origin  of  the  p  plane,  yields 


r  p  (v-c)t 

m 

V  -  c 


y(x,t)  = 


1  p  .  I*  -jO 

'  m  VZ-c2 


0  <  x  <  ct 


ct  <  x  <  Vt  }  0  <  c  <  V  (2.7) 


Vt  <  x 


J 


f  P  (c  -  V)t 

1  ~  '  a  t  ::2 


m 


c  -  V 


c-  V 


0  <  x  <  Vt 


y (x,  t)  =  <  ~  •  Vt  <  x  <  ct  >  0  <  V  <  c  (2.8) 


Ct  <  X 


From  (2.  7)  and  (2.  8)  the  velocities  are  readily  found  to  be 


z  .  S^Sl 

m  ,r2  2 


»  m 


V  -c 


V 

Tf - 1 


V  -  c 


0  <  x  <  ct 


ct  <  x  <  Vt  )  0  <  c  <  V  (2.9) 


Vt  <  x 


\ 

J 


8 


If  ■ 


JP 

m 


_P 

m 


(c  -V) 

~2 - 2 

c  V 


~T-~1 
c  -  V 


0  <  x  <  Vt 


Vt  <  x  <  ct  >  0  <  V  <  c  (2.  10) 


ct  <  x 


An  isolated  string  element  of  length  Ax  with  a  force  P 
traveling  over  it  in  a  time  At  receives  an  impulse  PAt  .  Let  the 
velocity  acquired  be  v  .  Then  the  momentum  acquired  is  m Axv  = 
PAt.  Since  At  =  Ax/V  ,  the  velocity  of  the  isolated  element  is 
v  =  P/mV.  Hence  the  impulse  delivered  is  P/Y  .  If  P  increases 
proportionately  with  V,  i.e.,  P/V  =  f  where  I  is  a  constant,  a 
series  of  these  elements  will  approximate  a  finite  length  of  string 
subjected  to  an  impulse  I  per  unit  length  which  imparts  a  velocity 
v  .  One  requirement  for  good  simulation  is  therefore  that  y Jv 
should  be  approximately  unit  (y^  3  &y/dt). 

Figures  2.  3(a)  and  2.  3(b)  show  the  velocity  and  displacement 
distributions  along  the  string  according  to  (2.9)  and  (2.  7)  with 
V  a  nc  (n  >  1).  The  force  ie  moving  supersonically  relative  to  the 
string  wave  velocity.  To  illustrate  the  degree  of  simulation  consider 
the  value  n  *  19.  From  Fig.  2.3(a)  it  ia  ssen  that  the  normal 
velocity  of  the  string  over  one  -nineteenth  of  the  distance  't  averted 
by  the  load  ie  5  percent  less  than  that  due  to  the  same  impulse  uni* 
formly  distributed.  Over  the  remaining  distance  it  is  about  0.  3 
percent  greater. 

Good  simulation  also  requires  small  displacements.  Figure 

2.3(b)  shows  the  displacements  with  a  maximum  of  y  =  vtn/(n+l)  or 
y/x  =  (v/V)n/(n+l)  *  y^/V  where  x  is  the  load  position.  Hence  the 
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V*ne 

n>l 


n*-l 
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(Q)  VELOCITY  DISTRIBUTION  —SUPERSONIC 


ratio  y  /V  or  v/V  should  be 
small.  In  order  to  estimate 
reasonable  values  of  the  ratio 
v/V  one  can  equate  the  kinetic 
energy  imparted  to  the  plastic 
work  done,  assuming  this  to  be 
much  larger  than  the  elastic 
strain-energy  capacity.  If  the 
final  strain  is  e  ,  the  yield  stress 


1 

- , 

i 

i 

C*flV 

n>l 

n 

J_ 

“1*. 

n+l 

♦p 

0 

l/n 

1  m  et 

is  0  ,  and  the  cross-sectional 

y 

area  is  A  ,  the  energy  equation 
is  mvCIZ  s  a  eA  .  Now  m  =  pA  , 

y  1/2 

(b)  DISPLACEMENT  DISTRIBUTION-SUPERSONIC  and  cy  =  (Oy/p)  is  the  maxi¬ 
mum  wave  velocity  of  the  string, 
p  being  the  mass  density.  Thus 
v  =  and  v/V  =j2e  /gy  where 

0y  -  V/cy  .  As  a  practical  example 
consider  an  aluminum  3tring 
stretched  almost  to  yielding  at 

o  -  50,  000  lb/in^.  With  a  mass 
Y  2  4 

density  p  =  0.  00025  lb-sec  /in 

the  wave  velocity  is  about  cy  = 

0.  355  mm/psec.  As  an  example 

of  an  explosive  with  one  of  the 

slower  detonation  velocities,  oxy- 

acetylene  gas  (50/50  mixture  by 

volume)  has  V  =  3  mm/psec  so 

that  0  as  9  .  Taking  a  large  strain 

Hence  the  maximum  displacement 


(cl  VELOCITY  DISTRIBUTION- SUUSONIC 


e«nV 

«>l 


o'  t/n  I 

Id)  DISPLACEMENT  DISTRIBUTION -SUBSONIC 

M-4WW-* 

FIG.  2.3  VELOCITIES  AND  DISPLACEMENTS 
OF  STRING 


of  e  =  0.08  the  ratio  v/V  =  0.045. 


is  approximately  y  =  0.  045  x  ,  where  x  is  the  load  position,  or  about 
5  percent  of  the  distance  traveled  by  the  load.  Using  the  sheet  explosive 
F,L  5Q6D  with  a  detonation  velocity  of  7  mm/psec  this  relation  becomes 
y  =  0.  019x  and  8y  20  . 
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Figures  2.3(c)  and  2.  3(d)  show  the  velocity  and  displacement 
distributions  along  the  string  according  to  (2.  10)  and  (2.8)  with 
c  =  nV  (n  >  1).  For  no  value  of  n  is  the  line  y^/v  =  1  approximated 
for  0  sx/ct  s  1/n  so  that  simulation  is  net  possible  when  the  load 
moves  subsonically  relative  to  the  string  wa^e  velocity. 

2.  3  MEMBRANE  THEORY 

An  infinite  stretched  membrane  is  subjected  to  a  ring  load  of 
magnitude  P  per  unit  length  of  circumference  (Fig.  2.2).  The  radius 
of  the  ring  increases  at  a  constant  velocity  V  .  Choosing  the  origin 
of  the  radial  coordinate  r  at  the  starting  or  detonation  point  when  the 
time  is  t  =  0  ,  the  equation  of  motion,  initial  conditions  and  boundary 
conditions  are 


6(r  -  Vt) 


(2.11) 


y(r.O)  =  (r,0)  =  0 


(2.  12) 


(0,  t)  =  0  Lim  y(r,  t)  =  0  (2,  13) 

r  -4  co 

where  y  is  the  deflection,  m  the  mass  per  unit  area  of  membrane 

1/2 

and  c  =  (S/m)  is  the  membrane  wave  velocity.  S  is  the  tension 
per  unit  edge  length. 

Let  y(X,t)  be  the  Hankel  transform  of  order  zero  of  the  function 
y(r,  t).  Then  applying  such  a  transformation  to  Eqs.  (2.  11)  and  (2.  12) 
gives,  with  the  aid  of  (2.  13) 

+  c2^y  ■  ™  J0^vt)  (2.  14) 

dt  mu 

y(x,°>  =  gf  (x*°) s  o  (2.  is) 


li 
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In  (2.  19)  P  |/2  an<*  Q  2/2  are  a8SOCiated  Legendre  functions 
of  the  first  and  second  kind  and  they  are  related  to  complete  elliptic 
integrals  of  the  first  kind  through  the  following  relations. 


P-1/2(A)  =  2F1  (1/2'  1/2:  =  2/ttK(/(1-A)/2  ) 

Q-l/2(  -A)  =  tt  /^/-2A  •  2F1(l/4.  3/4;  1;  1/A2)  = 
,/2/<l'-A)  K(/2/('l-A)  ) 


-1  <  A  <  1 


(2.21) 


-*  <  A  <  - 1 


where  is  a  hypergeometric  function. 


Before  applying  the  result  (2.  19)  to  evaluate  (2.  18)  the  latter 
should  be  rewritten  in  the  form 


z(a, 


J  Jo(ax)  Jo(8px)sinx  (l-p)dx 
0 


(2.22) 


In  the  first  double  integral  of  (2.22)  0  *  Bp  st  a,  while  in  the  second, 
a  £  8p  s  0  •  Making  use  of  the  relation  (2.21),  the  result  (2. 19) 
becomes 


f 0<l-p<8p-a  0<a<Bp 
\0<i-|i<a-8»i  0<8p<a 


/•  J  1  vrie  s  f0p-a<l-»*<8H+a  0<a<8p 

J  Jo(ax)Jo(0px)sin(i-p)xdx  a  <  TT(a8K)l/2  '  l'  \a-8p<i-HKa+S>4  0<8n<a 


■  ■  ■  ■>  ,«  k-K(k?)  a  +  <  1  < ' 


(2.23) 
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where  the  moduli  of  the  complete  elliptic  integrals  K  are 


1/2 


and 


in  which 


kj  =  l(l-A)/2] 
k,  =  [2{1-A)]1/2 


A  =  [a2  +  82  (a2  -  (l-|i)4]/2a&p 


from  (2.  20). 

Substituting  the  appropriate  result  from  (2.23}  in  the  integrals 
(2.  22)  leads  to  the  following  results 


1  <  a  <  0  TT(O0)1/2  a  =  Ij  (|i{.)  -  Ij  (^-) 


0  <  a<  1  <  0  TT/nox1/2-  -  T  tl+0-\  T  / 1  “tt \  i  T  /I -a . 
0  <  a<0  <  1  n<ae)  z  ’  ll  {T+g')  ’  h  (W)  h  {W} 

0  <  B  <a  <  I  tW/2.  *  I,  )  -  I,  +  I2  (^| ) 


(2.24) 


where 


t* 

IjGO  =  J  H  K(kj)dn 


I2(H)  =  J  H  k2K(k2)dp 
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It  is  of  interest  to  give  the  velocity  expressions  for  the  center 
of  the  membrane  since  they  can  be  obtained  explicitly.  Setting  r  =  0 
in  (2.  17)  and  continuing  the  derivations  as  before  leads  to 


1 1  -{tt/2  -sin'1  (1/8)}  /(@2-l)1/2]/(82-l) 


-  < 


1/3 


[  {cosh'1  (1/0)}  /( 1-82)1/2  -  1] /(I  -  a2) 


e  >  i 
8  =  1 
8  <  1 


(2.27) 


Figures  2.4(a)  and  2.  4(b)  show  the  velocity  displacement 
distributions  over  the  membrane  according  to  (2.24)  and  (2.  26)  when  the 
force  is  moving  supersonically  relative  to  the  membrane  wave  velocity. 
The  diagrams  have  been  drawn  for  the  case  V  =  2c  (8=2)  but  the  forms 
are  similar  for  all  the  supersonic  cases  (8>  1);  the  higher  the  value  of 
6  ,  the  flatter  the  velocity  curve.  Figure  2.  5  shows  the  velocity  for 
8=5.  As  in  the  representation  of  velocity  and  displacement  distri¬ 
butions  for  the  string  (F:g.  2.3)  the  velocity  v  used  to  render  results 
dimensionless  is  the  velocity  that  would  be  acquired  by  all  the  elements 
if  they  were  disconnected  from  each  other,  that  is,  v  =  P/mV  .  Hence 
good  simulation  of  impulse  applied  instantaneously  over  the  circular 
area  swept  out  by  the  detonation  front  is  obtained  if  the  curve  y(/v 
approximates  the  line  y(/v  =  1  (y^ady/dt)  .  Additionally,  for  good 
simulation  the  displacements  acquired  while  the  load  is  acting  should 
be  very  small. 


Figures  2. 4(c)  and  2.4(d)  shew  the  velocity  and  displacement 
distributions  over  the  membrane  according  to  (2.24)  and  (2.27)  when 
the  force  is  moving  subsonically  relative  to  the  membrane  wave 
velocity.  The  diagrams  have  been  drawn  for  the  case  V  ■  c/2  (8=  1/2) 
but  the  forms  are  similar  for  all  the  subsonic  cases  (8<  1).  Because 
of  the  significant  disturbances  running  ahead  of  the  toad,  good  simula¬ 
tion  is  never  possible  with  subsonic  loads.  Figure  2.4(c)  shows  that 
the  line  y  /v  -  1  is  not  approximated. 
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(b)  DISPLACEMENT  DKTWBUTCN -SUPERSONIC 


o  oa  w  to 

(0)  DISPLACEMENT  DWTR«UTWN-SUWONC 

m  mm  m 


FIG,  2,4  VELOCITIES  AND  DISPLACEMENTS 
OF  MEMBRANE 

is  then  about  14,  000  in/sec  or  about 


In  all  supersonic  cases  the 
maximum  velocity  occurs  under 
the  load  at  r  =  Vt  and  the  minimum 
velocity  is  the  velocity  in  the  region 
0  <r<ct  .  The  maximum  and  mini¬ 
mum  values  of  the  ratio  y^/v  are 
respectively 

82/(S2-l) 

and 

[  l-{n/2-sin‘1(l/8)}/{e2-l)1/2] 

b2/(b2-i) 

Table  2. 1  lists  values  of 
y(/v  for  several  values  of  0  and 
shows  the  extent  of  the  simulation. 
Although  the  minimum  values  of 
yt/v  require  high  values  of  B 
before  they  approximate  unity,  the 
radius  of  the  central  portion  of 
membrane  moving  at  this  velocity 
is  1/B  times  the  radius  of  the 
loading  circle. 

As  a  practical  example 
consider  an  aluminum  membrane 
with  a  mass  density  p*  0.  00025 
lb  -sec2 /ip4  stretched  to  a  stress 
of  o  *  50,000  lb/in2.  The  wave 
velocity  c  *  (S/m)*^2  *  (c/p)1^2 
0.  355  mm/psec.  The  detonation 


velocities  of  oxyacetylene  gas  (50/50  mixture  by  volume)  and  sheet 


explosive  EL  506D  are  approximately  3  mm/psec  and  7  mm/psec 
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a*r/ct 


FIG.  2.5  VELOCITY  DISTRIBUTION  ON  MEMBRANE  /9  -  V/c  -  5 


giving  values  for  Q  -  V/c  of  about  9  and  20  respectively.  From 
Table  2.  1  the  minimum  values  of  y(/v  corresponding  to  3=10  and 
20  are  low  by  about  14  percent  and  7  percent,  but  exist  only  in  central 
circles  of  radii  l/9th  and  l/20th  of  the  loading  circle.  The  less  the 
initial  stretching  of  a  membrane  (ox  string)  the  higher  is  the  value  of 
3  ,  for  a  given  explosive,  and  hence  the  better  is  simulation. 


To  estimate  the  deflections  which  may  exist  while  the  load  is 
still  acting  on  the  structure  the  initial  kinetic  energy  is  equated  to  the 
final  plastic  work,  assuming  the  latter  much  greater  than  the  elastic 
strain -energy  capacity  of  the  membrane.  If  each  element  has  a  final 


symmetric  strain  of  e  ,  the  yield  stress  is  0  and  the  membrane 

y  2 

depth  is  d  ,  then  the  energy  equation  is  approximately  mv  /2  =  2o  cd. 

1/2  y 

Now  m  *  pd  ,  and  cy  =  (oy/  p )  is  the  maximum  wave  velocity  in 
the  membrane,  so  v  =  CyV*5e  and  v/V  =v^c”/@y  where  0  =  V/cy  . 

For  a  strain  of  4  percent  and  By  =  9  the  ratio  v/V  *  0.  045  so  that, 
approximately,  the  deflection  is  y  =  0.  045  r  where  r  is  the  radius 


of  the  loading  circle.  The  larger  the  value  of  0^  the  smaller  the 
initial  deflections  and  hen  s  the  better  the  simulation. 
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Table  2.  1 

MEMBRANE— MAXIMUM  AND  MINIMUM  VALUES  OF  y %h 


0  =  V/c 

0  >  1 

(Vv)ma* 

=  e2/(e2-i) 

a  =  0 

Vv,min 

0  <  a  <  l 

2 

1.  233 

0.  527 

5 

1.  042 

0.751 

10 

1.  010 

0.861 

20 

1.  003 

0.  926 

50 

1.000 

0.969 

100 

It 

0.985 

200 

It 

0.992 

500 

II 

0.  997 

1000 

II 

0.998 

2. 4  CONCLUSIONS 

It  can  be  concluded  from  the  results  of  the  above  analyses  for 
strings  and  membranes  that,  for  good  simulation  of  distributed 
impulses  with  explosives,  the  detonation  velocities  V  must  be  greater 
than  the  wave  velocities  c  .  The  higher  the  velocity  ratio  V/c  ,  the 
better  is  the  simulation.  For  both  the  string  and  membrane  values  of 
V/c  greater  than  20  give  very  uniform  "initial"  velocity  distributions 
and,  provided  the  required  final  plastic  strains  are  not  too  large,  gives 
small  "initial"  displacements.  Values  of  V/c  greater  than  20  are 
certainly  quite  practical. 

The  actual  uniformity  of  the  initial  velocity  distribution  can  be 
seen  in  Fig.  2*  3(a)  for  the  string  and  in  Table  2.  1  (with  aid  of 
Figs.  2.4(a)  and  2.  5)  for  the  membrane.  Displacements  of  the  string 
and  membrane  acquired  during  loading  may  be  estimated  by  the 
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formulas  y  -  2ex/0 ^  and  y  =  4er/g^  respectively  where  e  is  the 
final  strain,  x  and  r  are  the  load  positions,  and  8^  =  V/cy  is  the 
maximum  wave  velocity. 

The  results  in  (2.  l]  and  [2.2]  indicate  that  it  is  desirable  to 
detonate  explosives  away  from  supports  to  minimize  their  influence 
while  the  load  is  being  applied. 
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3.  DYNAMIC  BUCKLING  OF  RECTANGULAR  PLATES 
IN  SUSTAINED  PLASTIC  COMPRESSIVE  FLOW 


3.  1  INTRODUCTION 


Dynamic  buckling  during  sustained  plastic  flow  has  been  described 
previously  for  cylindrical  shells  [  3.  l]  under  inward  radial  impulse  and 
for  rods  under  axial  impact  [3.2].  In  this  section  an  analysis  is  pre¬ 
sented  for  rectangular  plates  under  sustained  in-plane  flow,  and  the 
predicted  wavelength  is  shown  to  be  in  reasonable  agreement  with 
experimental  results. 

3.2  PRINCIPAL  FEATURES  OF  THE  PROBLEM 


A  simple  rectangular  plate,  flat  or  nearly  so,  is  made  to  flow- 
plastic  ally  for  a  short  time  by  sufficiently  large,  uniform,  in-plane 


stress,  with  three. components  0 


11  ’ 


°ZZ  ’  °12 


(Fig.  3.  1).  There  are 


non-zero  plastic  strain  rates,  related  to  the  stress  through  the  Prandtl- 
Reuss  equations  ([3.3],  p  39) 


e 


P 

11 


12 


(3.1) 


elastic  deformations  being  ignored.  Primes  denote  deviatoric  com¬ 
ponents  throughout. 


FIG.  3.1  NOTATION 


While  this  in-plane  flow  is  proceeding,  there  may  also  be  small 
flexural  motion  in  the  transverse  (x^)  direction,  resulting  in  a  plastically 
buckled  form.  This  flexure  alters  the  strain  rates,  except  on  the  middle 
surface.  But  it  is  postulated  that  it  does  not  cause  reversal  of  stiaiii* 
rate,  so  we  are  concerned  with  loading  only,  in  the  sense  of  plasticity 
theory.  This  is  the  type  of  flexural  buckling  envisaged  by  Shanley  [  3.  4] 
for  the  column,  and  considered  for  the  plate  under  uniaxial  compression 
by  Pearson  [3.5]. 

In  particular  the  transverse  motion  of  the  plate  introduces 
differences  between  the  strains,  and  strain  rates,  of  a  point  on  the 
upper  face  of  the  plate  and  the  point  below  it  on  the  lower  face.  Then 
the  strain-increment  vector  of  the  flow  rule  is  not  the  same  for  the  two 
points,  and  in  general  the  directions  will  be  different.  The  representa¬ 
tive  points  on  the  appropriate  current  yield  surfaces  must  then  be  such 
that  the  normals  are  correspondingly  different,  implying  a  stress 
difference  which  can  appear  even  if  strain  hardening  is  absent.  This 
can  not  occur  for  the  one -dimensional  stress  of  the  column.  If  there 
is  no  strain  hardening,  flexure  of  the  column  can  occu.'  without  inducing 
bending  moment.  In  the  plate  the  stress -differences  imply  bending 
moments  and  also,  in  general,  twisting  moments. 


The  strain  rate  for  <i  lower -face  point  maybe  expressed  by 

it?.  +  Ac?.  ,  and  for  the  upper -face  point  by  e  Pi  -  Ac?.  .  The  plastic 
iJ  iJ  .  ij 

parameter  is  X  for  the  middle  surface,  but  X  +  AX  for  the  lower-face 
•  ■ 

point  and  X  -  AX  for  the  upper-face  point.  With  a  similar  notation  for 
stress,  the  Prandtl-Reuss  equations  for  the  lower  face  are 


e?  +  Ac?  = 


{X  +  AX)(o.'.  +  AOj'.  ) 


and  for  the  upper  face  are 


(3.2) 


(X  -  AX  ){o'  -  Ac  '  ) 

J  *  J 
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We  take  these  for  the  in -plane  directions  only  (i,  j  =  1,2).  Eliminating 
e?.  we  have 

=  (X  +  AX)(o.'.  J  AO.'.)  -  (X  -  AX) (a/.  -  AO. ' ) 

The  difference  quantities  however  are  to  be  small,  corresponding  to 
small  flexure,  which  must  prevail  at  least  in  the  earlier  stage  of  the 
motion  of  a  nearly  flat  plate.  These  relations  will  accordingly  be 
linearized  to 


«fj  =  + 

consist  of 

Ae^j  =  XAOjj  +  AX 

(3.3) 

A®22  =  °22  +  °ZZ 

(3.4) 

Aej2  =  ^Ao,2  +  a12  AX 

(3.5) 

By  elimination  of  AX  ,  (3.  3)  and  (3.  4)  yield 


Ae 


:T? 


11 


Ae 


A', 


11 


45  - 

a22  an 


^°22 

°Zl 


(3.6) 


Adding  (3,3)  and  (3.4),  we  obtain  a  symmetric  equatic.'.  Combining 
this  with  (3.5)  to  eliminate  AX  ,  we  have 


*‘11+^22 
°\\  +  022 


Ae 


12 


=  X  ( 


12 


A  Jl  1  +  Ag22 

°n  +  a22 


A  a 
a 


12 


(3.  7) 
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The  unperturbed  flow  is  regarded  as  given.  This  means  that  the 
strain  rates  in  (3.  1}  are  given,  as  well  as  the  stress,  and  therefore 
\  is  given.  It  will  be  evident  that  the  effects  of  strain  rates  introduced 
by  some  buckling  deflection  w(Xj,  x^,  t)  on  strain-increment  vector 
directions  will  differ  according  as  the  original  strain  rate  is  large  or 
small .  The  original  strain  rate  is  therefore  a  necessary  part  of  the 
specification  of  the  state  from  which  buckling  occurs. 


Now  (3.  6)  and  (3.  7)  provide  two  linear  equations  relating  the 
unknown  stress  differences  AOj  j  .  Ao^  •  ^a\z  *  to  *train-rate 
differences  Ae^j  ,  A j2  • 

In  the  next  section  we  obtain  a  third  linear  equation  from  the  strain¬ 
hardening  law.  The  three  equations  will  then  play  a  part  similar  to  that 
of  moment -curvature  relations  in  elastic  plate  theory.  Combination  with 
the  equations  of  motion  of  the  plate  element  then  results  in  a  differential 
equation  for  the  deflection  w  alone. 


Experimental  results  reported  later  in  this  paper,  and  interpreted 
by  means  of  the  present  theory,  ar .  all  such  that  the  observable  defor¬ 
mation  is  predominantly  plastic.  For  this  reason,  and  to  minimize 
complications,  we  now  assume  that  elastic  deformation  can  be  entirely 
neglected.  Then  the  superscript  p  can  be  dropped  in  (3.  1)  through 
(3.  7),  and  we  can  further  write 


Ac 


11 


h 

7 


0  w 

~T2 

3x, 


Ae 


22 


v,  a2* 
h  o  w 

?1*2 


Ar 


12 


h 

133? 


A2* 
o  w 


0.8) 


If  the  plate  is  not  initially  fiat,  the  initial  ordinate  of  the  mid* 
surface  will  be  denoted  by  »fJx,)  ,  and  then  w  means  the  observable 
deflection,  the  additional  ordinate,  at  time  t  . 
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3.  3  ADAPTATION  OF  A  ST  RAIN -HARDENING  LAW 


The  strain-hardening  law  to  be  employed  here  is 


where 


a 


I 


H  (  de) 


(3.9) 


de  ^/|(de..  de..)  ;  i,  j  =  1,2,3  (3. 


10) 


define  the  positive  "equivalent  stress"  and  "equivalent  strain  increment" 
corresponding  to  [  3.  3]  p  26  Eq.  (14),*  and  p  30  Eq.  (16),  except  that 
the  present  notation  omits  bars.  The  integral  in  (3.8)  is  path -dependent, 
and  the  strain  paths  are  not  the  same  for  different  elements  of  material 
on  the  same  thickness  line. 


In  the  present  problem  we  regard  the  transverse  shear  strains 
C23  *  e31  38  negligibly  small,  supposing  that  thickness  lines  remain 
normal  to  the  middle  surface  and  straight,  as  commonly  done  in  thin 
plate  and  shell  theory.  Then,  in  view  of  plastic  incompressibility, 
we  can  write  from  the  second  of  (3.  10), 

(dc)2  =  yl(deu)2  +  (de2)2  +  dCj  .  de2  +  (deu)2]  (3.  11) 

The  differentials  in  (3.  10)  are  now  taken  specifically  as  those  pertinent 
to  the  midaurface,  and  occurring  in  a  time  increment  dt  . 

For  the  lower  surface  the  corresponding  strain  increments  are 

d,„-!siiit  •  •  dti2-riud*  <312> 

where 

Xu  a  32w/8x2  ,  i22  =  d2w/dx2  ,  r.J2  =  d2w/dxjdx2  (3.13) 
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Writing  de,  for  the  equivalent  atrain  increment  at  the  lower 
l  2 

surface,  we  may  evaluate  (de^)  from  (3.  11)  by  introducing  the 
replacements  (3. 12). 

We  have  already,  in  deriving  (3.3),  (3.4)  and  (3.5),  treated  the 
difference  terms  introduced  by  the  flexure  as  small  in  comparison  with 
the  terms  of  the  dominating  flow.  Correspondingly  in  (3.  12)  the  flexural 
w  terms  will  be  supposed  small  compared  with  the  midsurface  terms. 
Then  the  squares  and  products  from  (3.  12)  can  be  linearized  with  respect 
to  w  .  We  find 

(d€j)  =  (de)  ’  jhdttdCjjfHjj +2-*22)  f  d€22^*22 +F*11^  +  de12*12^ 

(3.  14) 

2 

For  the  upper  surface  (de^)  is  obtained  by  changing  the  -  to  +  . 

We  apply  the  law  (3.  9)  separately  to  the  lower  and  the  upper 
surfaces,  and  wish  to  obtain  the  difference 


ZtC 


=/dct-/dcu 


(3.  15) 


which  defines  Ac  .  From  (3, 14)  and  the  companion  expression  for 
2 

(dc  )  ,  we  have 


But 


(dc^)2  *  (dcu)2  =  (dct  +  dcu)(dct-deu) 


(3.  16) 


=  -yhdttdCjjfKjj  +»r  ^2)  +  ••  •  ] 


de,  +  dc  =  2de 
t  u 


(3.17) 
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Introducing  this  in  (3.  16),  and  converting  to  rates  by  inserting  dt  as 
required,  we  find 


def,~deu  “  "Thdt  *  *  ^il^ll"^  *22^  +  ^22^*22 +?'H11^  +  ^12*12^  (3*  18) 


and  now  (3.  15)  can  be  written  as 
t 


Ae  =  "  ThJ^J^ll(^ll+2  ^22>  +  ®22^*22  +  2*1 1^  +  ^  12^12^  dt  ^3‘  ^ 
o 

This  becomes  much  simpler  if  the  ratios  e ^  ^ / e  ,  2 /s  »  e  ^/e  ,  per¬ 

taining  to  the  midsurface  flow,  are  independent  of  t  throughout  the 
motion.  In  general  the  Levy-Mises  relations, to  which  (3.  1)  now  reduce, 
imply 

e..  3c!. 

=  To  (3.20) 


We  now  suppose  that  the  imposed  stress  components  Ojj  ,  o ^ 
are  maintained  in  constant  ratio  throughout  the  motion.  This  implies 
constancy  of  the  right-hand  side  of  (3.20).  Then  the  ratios  Ejj/e  etc. 
in  (3.  19)  are  in  fact  constants,  and  the  integration  can  be  effected  by 
removing  the  dots  on  the  k's  .  Thus  (3.  19)  becomes 


Ae  =  -  h  [a(K11+2‘H22)  +  ®(K22  +  7*11^  +  Yh12  1 


(3.21) 


where 


2e 


a  = 


11 


3e 


°H  .  _  2c22 

“5  ■  *  • 

3e 


ll 


2c 


Y  = 


12 

3e 


12 


H11  =  a2w^xi 


etc. 


(3.22) 


and 


°U  a22 


a 


(3.23) 


The  quantity  2Ae  ,  expressed  by  (3.21)  in  terms  of  w  ,  shows 
how  much  further  along  the  equivalent-stress— equivalent-strain  curve 
(3.  9)  the  point  for  the  lower  surface  is  in  advance  of  the  point  for  the 
upper  surface.  The  difference  in  the  values  of  0  for  the  two  surfaces 
is  given  by  the  differential  form  of  (3.  9), 

2A0  =  2H '  Ac  (3.24) 

and  we  shal.  take  H  '  to  be  constant. 

Again  with  reference  to  the  experimental  results  we  seek  to 
explain,  and  for  brevity,  we  now  restrict  the  analysis  to  the  case 
°12  =  ®  *  Then  the  first  of  (3.  10)  yields 

2o  Ao  =  (2on  -  o22Uan  +  (2 o22  -  °u)Ao22  (3.25) 

Now  (3.  24)  becomes,  with  (3.25)  and  (3.21) 

(2air°22Uoll  +  ^2o22  "°ll^a22  =  "2H  oh^  a(Hli  +2"k22^ 

(3.26) 

+  0  (h22  +  2" K  i  i  ^ 

This,  and  (3.  6)  (3.  7),  form  a  set  of  three  equations,  which  are  linear 
relations  between  the  stress  differences  AOjj  ,  Ao22  ,  AOj2  ,  and  the 
derivatives  of  w  [in  view  of  (3.8)]. 

Since  we  contemplate  small  stress  differences  occupying  only  a 
small  segment  of  the  line  represented  by  (3.  9),  and  there  is  no  un¬ 
loading,  the  variation  of  stress  through  the  thickness  may  be  taken  as 
linear.  Then  the  moments  per  unit  run  ,  M22  ,  M^2  are  given 

by 

h2A0u  =  6Mjj  ,  h2A  o22  =  6M22  ,  h2ACj2  =  6Mj2  (3.27) 
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We  rewrite  (3.  6),  (3.7)  and  (3.26)  using  these,  and  (3,8). 
(3.  6)  becomes 


M 


12  =  “ 


Then 


(3.28) 


(3.  7)  becomes 


°22M11  '  allM22 


12X 


T  [  <2o22  “  CT1 1)  H1 1  "(2all  ‘°22)k22^  (3*29) 


and  (3.26)  becomes 


(2al l”°22^Mll  +  ^°ZZ~  all^M22  =  “TH  0h  ^Kll  +  TrzZ^  +^K22+<rKU^ 


(3.30) 


In  (3.  30)  we  can  put  H#  =  0  for  a  perfectly  plastic  material,  and  still 
obtain  non-zero  Mjj,  from  (3. 29)  and  (3.  30).  These  represent 

the  stress  differences  induced  by  movement  of  the  stress  points  on  the 
now  unchanging  yield  surface,  referred  to  in  Section  3.  1,  and  depending 
on  the  curvature  rates. 


The  solutions  of  (3.29),  (3.30),  for  Mjj  ,  •  *re 


~Mn 

.3 

-  .  f  (In  n  V>n  -a  1 

poJ2-ouy 

M22 

7  1  '*°22  CT1 1  ,K1 1  '  ll  aZZ,KZZl 

■ 

L_  -J 

‘TT"  [a(Kll  +rK22)  +  S{*22  +  Z)] 


J11 


3ZZ 


(3.31) 


We  shall  xefer  to  the  terms  of  the  1st  line,  having  the  factor  1  / X  , 
as  the  "directional”  moments,  and  to  the  terms  of  the  2nd  line,  pro¬ 
portional  to  H  / ,  as  the  strain-hardening  moments. 
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3.4  EQUATIONS  OF  MOTION 


A  plate  element  dx  dy  is  subjected  to  uniform  stress  ,  c ^  , 
together  with  the  moments  Mj(  .  M22  ,  Mj2  »  and  transverse  shear 
forces  per  unit  run  Qj  ,  necessarily  associated  with  these.  We 
neglect  rotary  inertia.  The  equations  of  motion  take  the  well-known 
forms  of  plate  theory 


dQi  2Q2 
33T  *  37  +  °llh 


*2  J- 

72  *  °22n  rr  vh  r  z 

3  x  j  ilx2  3 1 


_  *Mn  SM12 
l  3x2 


3M22  .  dMl2 
2  =  + 


where  p  means  density. 

\ 

By  elimination  of  Q  ,  ,  Q-, 

*  M 

.  ,  a*M12  .  S  U2?.  . 

'7~T  +  2777T,  +  77— r  °nh 

Ojtj  12  dx. 


32(w+w_)  32(w+w  ) 

— *  °22h  — rT~ 

®*2  3*2 


ii.32) 


In  this  w*  now  use  (3,  31)  for  Mjj  .  M22  ,  icr  ^12  ’  convert 

it  to  a  differential  equation  for  w  alone  to  obtain 


.2  f  2o,>-olt,  (2c. o,,H2o,^-o, .) 

h  J,  22  11*2  v?  w  A  ,  il  2242  3  w  U  22'  22  11' 

- T  \  i  — r— —  >  — ~  ~r  *  l  “V ' —  f  7  *  t  *  11  2  'r" 

24X  1  0  3 ibi*  u  3t3xJ  <T 


*  j  °ll  ,  1  flV  34w  °22  ,a  .  I  .  34w 

t~~  y-r(a4rBi~4 *  —  (8+?a,^r 


1  ‘-(S  ♦  ta )  ♦  ^  (a  *  ))  •  C, 


4-  l  ♦  ya )  4  (a  *  j-S  )]  -£ 


^(*4^1  *>♦*  ) 

l  ~T“*  *  °22  — ~'T~L 
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1  2 j 


•  cr 
3 1 


(3.33) 
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All  the  coefficients  are  constants,  except  for  X  in  the  first 
line.  As  we  have  observed  in  Section  3.  1,  X  must  be  specified,  and 
this  is  equivalent  to  t,?ecifying  imposed  strain  rates  Cjj  ,  as 

well  as  imposed  stresses  Cjj  ,  •  Since  we  have  postulated  slight 

strain  hardening,  it  is  reasonable  to  treat  0^  ,  c  ^  as  constant 
throughout  the  flow.  We  now  postulate  farther  that  the  imposed  strain 
rates  are  to  be  constant,  and  therefore  X  will  also  be  constant.  Now 
(3.  32)  is  a  linear  equation  with  constant  coefficients.  The  uolution  for 
a  rectangular  plate  with  simply  supported  edges  is  examined  in 
Section  3.  5  below. 

Since  from  (3.  22) 


11 


=  2a  +  6 


'T 

22 


a  +  26 


wc  v  an  rewrite  {3.  33)  as 
2 


_h 

24 


X  Zt  dx? 


ZairJ22  2  d5w 


dt  dx 


[“2<2oU'C22,(Z02?.*°11 
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3^ 


+•  o* 


22 


(3.  34) 


ox. 


If  the  compressive  flow  and  the  buckling  are  occurring  slowly,  the 
inertia  term  may  be  negligible.  Eut  it  is  not  evident  that  the  first  line 
of  terms  may  be  dropped,  since  X  appears  in  the  denominator  and  is 
small  tor  a  small  flow  rate.  The  magnitudes  of  the  time  derivatives 
of  w  will  evidently  depend  on  the  initial  displacement  and  velocity. 
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3.5  THE  SIMPLY  SUPPORTED  RECTANGULAR  PLATE 
We  take  o 


'll 


a  f  a  ,,  =  0  ,  corresponding  to  uniaxial 

°  2  "2 


comprespive  flow.  Then  o  =  .  The  differential  equation  (3.34) 

reduces  to 


v2  a5 

h  o  w 


*5 
o  w 


.5 
0  w 


- 4  t 

24X  dtSxj 
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atdx2 

o - y 

St 

77 
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*2 
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°o7Z  +  p 
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<v 

X 

h*  ro 

(3e35) 


If 


w  =  I  Ea  sin  — 

o  mn  a 


mTTX  HTTX- 

1  sin 


(3.36) 


we  may  take 


mnx .  ntTX? 

w  =  II  W  (t)sin - -  sin  — r — 

mn'  a  b 


(3.37) 


This  makes  in  (3.  31)  zero  at  the  edges  Xj  *  0,  a;  x2  =  0,  b, 

so  the  edges  are  simply  supported.  Then  (3.35)  requires 


K  lisn.)4  +  4(^i4  +  b  (EL>2  if )2  ] 

12  X 


+  ( 


21)2[ 


1  hZH'(-22.)Z-8]W  +  OW 

&  o  mn  mn 


T2 


s  a  (S3)2 
o  mn'  a  ' 


3.38) 


For  suitably  limited  initial  deflection  wQ  and  velocity  w(Xj,  x.,  ,  0) 

the  time  derivatives  of  W  can  be  vary  small.  These  would  be  the 

mn 
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circumstances  of  a  "static"  test  with  0q  increasing  rather  slowly  to  a 
critical  value.  The  result 


W  = 
mn 


o  a 
o  mn 


>V!=H!)2 


0 

O 


(3.39) 


then  indicates  a  critical  compressive  stress 

°o  =  72  b2^^)2  (3.40) 

It  is  remarkable  that  this  is  independent  of  b/n  ,  the  half -wavelength 
in  the  direction  transverse  to  the  load.  For  m  =  1  it  is  the  simple 
tangent  modulus  formula. 


The  same  critical  condition  is  obtained  from  (3.  35)  if  we  take 
w  in  the  form 


mrrXj 

w  =  W  sin - 

m  a 


(3 


4  I  t 


i.e,,  independent  of  . 


Since  =  0  ,  the  value  of  given  by  (3.30)  consists  of  the 

directional  part  only,  and  is  therefore  very  small  in  the  "static"  test, 
for  (3.  37)  or  for  (3.41). 


Returning  to  the  dynamical  problem  of  ihe  simply  supported  plate 
with  the  deflection  (3.  36),  we  suppose  that  the  imposed  strain  rate 
Gjj  is  large  compared  with  any  additional  strain  rates  introduced  bv 
w  .  Then  the  directional  moments  in  (3.  31)  could  be  neglected  rom- 
pared  with  the  strain-hardening  moments.  The  first  line  of  (3.34),  and 
that  of  (3,  38),  disappear. 
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Considering  an  imposed  0q  much  larger  than  the  statical 
critical  value  (3.  39),  we  shall  have  (3.  38)  in  the  form 

W  -  -  ( J2HI  )2  [  j  -  -L  h2H'  (~)2]  W  =—  a  (-— ) 2  (3.42) 

mn  p  a  o  TZ  a  mn  p  mn  a  '  ' 


It  is  convenient  now  to  change  to  a  dimensionless  deflection  co¬ 
efficient  w  ,  a  dimensionless  time  t  ,  and  to  introduce  a  dimension- 
mn  ^ 

less  constant  s  ,  defined  by 
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w 


mn 


mn 


T  = 
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H 


TZp 


,  12o  a 

2  o 

8  =  ~2 — TZ 

tt  H  h 


(3.43) 


Then  (3.  42)  becomes 


,2 

d  w 


dT 


mn 

1 


2  .  2  2.  2  2  jnn 

m  (8  -m  )  w  =  s  m  -r— • 
mn  n 


(3.44) 


For  m  <  s  the  solution  will  contain  hyperbolic  functions  of  time 
which  may  become  large  during  the  motion. 

For  a  plate  initially  perfectly  flat  we  have  amn  =  0  ,  and  any 
deflection  is  due  to  initial  lateral  velocity.  Then  if 


mnx  j  nrrx^ 


w(Xj,  x2  ,  0)  =  0  ,  ^-wt(xi,x2,  0)  =  vo  EEemn8in~ — sin  —  (3.45) 


the  result  obtained  after  integrating  (3.  44)  is 


,  .  ®  r  i  mirx.  mrrx- 

fw(x  ,x  ,t)  =  V  (  E  E  —  S  sinh  p  T  •  sin— —  •  sin  -t- — 

h  1  2  o  )  _  t  i  p  mn  m  a  b 
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,  ,  p  mn  m 

n=  1  m=  1  m 


mux,  mrrx_ 

~  •  s  in  — ^  ( 3 .  46) 


where 


Pm  >  0  . 

m 


m^(s^-m^)|  ,  r  <s  <r  +  1  (3.4?) 

I 


m 
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This  shows  that  the  series  for  the  deflection  at  time  t  is  obtaineo  from 

the  series  in  (3,  45)  for  the  initial  velocity  by  applying  to  each  coefficient 

a  magnification  factor  (but  for  a  constant)  which  is  either  sinh  p^T 

or  pm*  sin  pmT  f,  the  former  being  potentially  the  greater.  Examining 

it  further  we  observe  that  if  s^  is  suitably  large  compared  with  unity, 

2 

p  as  a  function  of  m  will  be  small  for  m  =  1  and  for  m  close 
m 

to  s  (i.e.,  m  =  r),  but  will  have  a  maximum  at  m  -  s ,/yT ,  with  the 
value  Pm  max  =  s^/2  .  The  hyperbolic  magnification  factor  itself,  as 
a  function  of  m  ,  also  has  a  maximum  for  the  same  m  ,  and  this 
maximum  is 

(Pm‘  sinh  PmT  U  *  <3'48> 

s 

This  can  become  very  large  as  T  increases.  If  it  does,  there  is 
conspicuous  magnification  of  harmonics  having  m  in  the  neighborhood 
of  ,  regardless  of  n  . 

To  illustrate  the  magnitudes  involved,  we  consider  an  aluminum 
alloy  with 


o  =  2.  5  x  10^  psi  ,  H  '  =  1.  3  x  10^  psi 


(3.49) 


Then  from  (3.  43) 


2  1  a 

8  =  47T* 

h 


8  "  2j5T  IT  ’  m’  = 


1  V* 


dr  c  <3-50> 


Evidently  there  are  hyperbolic  terms  in  (3.46)  when  a/h  is 
greater  than  2.  07  (i.  e.  ,  s  >  1).  Much  larger  values  of  a/h  can  allow 
comspicuous  magnification  of  harmonics  if  (3.48)  can  become  much 
larger  than  unitv.  This  of  course  depends  on  the  duration,  the  range 
of  t  in  which  the  compressive  flow  stress  0q  is  maintained  and  no 
strain-rate  reversal  has  yet  occurred.  The  most  magnified  harmonic 
has  m  close  to  nij  as  given  in  (3.  50),  and  the  corresponding  half¬ 
wavelength  is 

—  =  2.94h  (3.51) 

ml 


37 


The  bar  problem  analogous  to  the  present  plate  problem  was 
treated  earlier  [ 3.2] .  The  most  magnified  wavelength  for  a  bar  of 
rectangular  section  of  thickness  h  is  the  same  as  (3.51).  Thus,  the 
hinges  on  the  boundary  conditions  do  not  affect  the  wavelength. 

If  the  initial  velocity  is  zero,  any  deflection  is  due  to  initial 

deflection  wq  .  Taking  wq  in  the  form  (3.36),  we  have  again  (3.42) 

for  W  in  (3.36),  or  equivalently  (3.44)  for  w  in  (3.43).  The 
mn  mn 

appropriate  solution  is  now 


mn 


w(Xj,x2,t)  =  a  Z  Z  ~2 

m=l  n=l  s  -m 


f  cosh  .  mr*l  .  rmx; 

2[cos  PmT-l3sm_  sxn 


(3.52) 


where  cosh  is  read  for  m  <  s  ,  cos  for  m  >  s  . 


The  magnification  factor  applicable  to  the  mn  term  of  wq  is 


now 


-2 - j  (ao sh  p^T  *  1) 

s  -  m 


(3.53) 


This  factor  has  occurred  previously  in  the  treatment  of  the 

bar  [3.2j.  Its  discussion  is  repeated  here  for  convenience.  Supposing 

that  s  is  so  large  that  there  are  several  hyperbolic  terms,  we  consider 

2 

(3.53)  as  a  function  of  m  .It  has  a  maximum  where 


s 

2  m' 


1  =  - 


cosh  PmT“  1 

p  t  sinh  p  t 
rm  m 


(3.54) 


Conspicuous  magnification  depends  on  cosh  p^T  becoming  much 


But  in  this  case  sinh  p  t  has  an  approximately 


larger  than  unity. 

equal  value.  Replacing  the  right-hand  side  of  (3.  54)  by  -1/p^T  ,  we 

may  then  observe  that  p  t  must  be  considerably  larger  than  unity, 

m  2 

so  1/p  r  is  small.  Consequently  (3.54)  will  require  a  value  of  m 
m  2  2 

slightly  greater  than  s  /2  .  T?king  s  / 2  as  an  approximation,  the 
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2 

greatest  magnification  fu (3.  53)  becomes  approximately  2  cosh  1/2  s 
As  was  shown  in  [3.2],  this  in  fact  becomes  large  in  durations  of  the 
order  of  lOOpsec. 

Since  the  magnification  factor  is  greatest  for  m  close  to  s/^/z  , 
the  most  magnified  wavelength  is  the  same  as  in  the  case  of  initial 
velocity  perturbations. 

3.  6  COMPARISON  WITH  EXPERIMENT 

Experiments  were  undertaken  in  which  square  tubes  of  aluminum 
(6063 -T5)  were  projected  end-on  against  a  massive  steel  plate.  Each 
side  of  the  tube  behaved  as  a  plate  in  uniaxial  compression.  To  facili¬ 
tate  acceleration  of  the  tubes,  they  were  mounted  on  round  rods  which 
could  be  fired  from  a  standard  rifle.  Several  specimens  are  shown  in 
Figs.  3.2  to  3.9.  As  indicated  in  the  figures,  several  combinations 
of  boundary  conditions  were  used. 

The  half -wavelengths  from ‘the  specimens  of  Figs.  3.2  to  3.9 
are  listed  in  Tables  3.  1  -  3.4.  The  1/ 16-inch-thick  plates  (Tables 


300  (SAC -2) 


400  i SAC  -  I  i 


OP  4WI  22 


FIG.  3.2  l/lfc-INCH  WALL  BY  5/8-INCH-SQUARE  ALUMINUM  TUBING  BUCKLED 
BY  AXIAL  IMPACT  WITH  A  MASSIVE  STEEL  PLATE 
AT  THE  INDICATED  VELOCITIES 
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3.  1  -  3.  3)  show  a  variation  in  half -wavelength  with  boundary  condition 
from  0.26  inch  for  1/2-inch-wide  plates,  hinged  or  clamped  on  the 
long  edges  (depending  on  the  response  of  the  adjoining  plates),  to  0.61 
inch  for  plates  which  are  free  on  the  long  edges.  As  indicated  by  the 
large  buckles  in  the  second  and  third  specimens  of  Fig.  3.2,  the  wave¬ 
length  does  not  appear  to  depend  significantly  on  whether  the  long  edge 
is  hinged  or  clamped.  The  1/ 8-inch -thick  plates  were  tested  with 
hinged  or  clamped  boundaries  only  on  the  long  edges  and  show  a  half¬ 
wavelength  of  about  0.45  inch.  From  (3.43)  the  theoretical  half¬ 
wavelength  X  corresponding  to  the  most  magnified  harmonic 
is 


184  f  t/ sec  ( i  AC  -  i ) 


344  U  AC  3) 


C.P  ?  3 


FIG.  3.3  1  /16-INCH  WALL  BY  7/8-iNCH-SQUARE  ALUMINUM  TUBING  BUCKLED 
BY  AXIAL  IMPACT  WITH  A  MASSIVE  STEEL  PLATE 
AT  THE  INDICATED  VELOCITIES 
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82  U/sec 


(2CSC- I ) 


(2CSC-2) 


GP  4999  24 


FIG.  3.4  1/16-INCH  WALL  BY  5/8-INCH-SQUARE  ALUMINUM  TUBING, 

WITH  TWO  OPPOSITE  CORNERS  CUT,  BUCKLED  BY  AXIAL  IMPACT 
WITH  A  MASSIVE  STEEL  PLATE  AT  THE  INDICATED  VELOCITIES 


1 1  8  ft/sec 


(  2  CSC  -  I  ) 


123  (2C5C-2)  or -499S 

FIG.  3.5  1/1 6" INCH  WALL  BY  7/8-INCH-SQUARE  ALUMINUM  TUBING, 

WITH  TWO  OPPOSITE  CORNERS  CUT,  BUCKLED  BY  AXIAL  IMPACT 
WITH  A  MASSIVE  STEEL  PLATE  AT  THE  INDICATED  VELOCITIES 


59  ft/sec 


(4C5C  -3) 


* 

GP’*?9<*  2  b 


FIG.  3.6  116-INCH  WALL  BY  5  8-INCH-SQUARE  ALUMINUM  TUBiNG, 

WITH  FOUR  CORNERS  CUT,  BUCKLED  BY  AXIAL  IMPACT 
WITH  A  MASSIVE  STEEL  PLATE  AT  THE  INDICATED  VELOCITIES 


0»  «**•*> 


FIG.  3.7  l  lo-INCH  WALL  BY  7  8-INCH-SQUARE  ALUMINUM  TUBING, 

WITH  FOUR  CORNERS  CUT,  BUCKLED  BY  AXIAL  IMPACT 
WITH  A  MASSIVE  STEEL  PLATE  AT  THE  INDICATED  VELOCITIES 


FIG.  3.8  1  8-  INCH  WALL  BY  3  4-INCH-SQUARE  ALUMINUM  TUBING  BUCKLED 
BY  AXIAL  IMPACT  WITH  A  MASSIVE  STEEL  PLATE 
AT  THE  INDICATED  VELOCITIES 


TRUE  STRESS 


As  indicated  in  Fig.  3.  10,  the  factor^/H  7oq  varies  considerably 
from  1  to  10  percent  strain,  decreasing  as  strain  increases.  Thus; 
(3.55  )  indicates  that  the  half -wavelength  should  decrease  as  strain 
increases.  This  expectation  is  borne  out  by  the  data,  as  may  be  seen 
by  comparing  Tables  3.  1  -  3.  3.  The  effect  of  freeing  the  long  edges 
appears  to  be  to  increase  the  instability  so  that  buckling  occurs  at 
lower  axial  strains. 


FIG.  3.10  STRESS-STRAIN  CURVE  FOR  TUBING  MATERIAL  (Aluminum  6063-T5) 


In  Table  3.  5  theoretical  half-wavelengths  are  given  for  values 
corresponding  to  strains  in  the  region  of  interest.  For 
the  1/  16-inch-thick  plates  the  predicted  half-wavelengths  for  strains 
of  1  and  5  percent  are  0.26  and  0.  14  inch.  The  observed  half¬ 
wavelengths  vary  from  0.26  to  0.61  inch.  The  difference  is  attributed 
to  the  effects  of  early  elastic  deformations  on  the  subsequent  plastic 
deformation  [3.4],  For  the  1/8 -inch -thick  plates  the  predicted  half¬ 
wavelengths  for  strains  of  1,  5,  and  10  percent  are  0.52,  0.27,  and 
0.  19  inch.  The  observed  half -wavelengths  (Table  3.4)  fall  within  the 
values  for  1  and  5  percent  strain. 

The  agreement  between  theoretical  and  experimental  wavelengths 
is  reasonable,  indicating  that  the  theory  describes  the  mechanics 
fairly  well.  The  largest  factor  yet  to  be  accounted  for  in  the  theory  is 
the  variation  in  hVo  with  strain. 


Table  3,  5 

CALCULATED  HALF-WAVELENGTHS 


h 

(inch) 

Assumed 

Buckling 

Strain 

<%) 

a 

0 

H ' 

X. 

K 

\ 

(inch) 

0.062 

1 

0.  096 

4.  13 

0.26 

5 

0.35 

2.  17 

0.  14 

0.  125 

1 

0.  096 

4.  13 

0.52 

5 

0.35 

2.  17 

0.27 

10 

0.69 

1.54 

0.  19 
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4.  DYNAMIC  BUCKLING  OF  A  THIN  CYLINDRICAL  SHELL 
UNDER  AXIAL  IMPACT 

4.  1  INTRODUCTION 

Experimental  investigations  by  Coppa  [4.  i]  and  others  have 
shown  that  thin  cylindrical  shells  subjected  to  axial  impact  can  buckle 
dynamically  into  the  familiar  diamond  pattern  of  large  deflection  static 
buckling.  The  diamonds  in  the  dynamically  buckled  shells,  however,  are 
smaller  than  in  the  static  pattern  because  higher  modes  respond  faster 
than  lower  ones  and  soon  predominate  the  motion.  Dynamic  buckling 
into  very  high-order  modes  (short  wavelength)  has  also  been  observed 
in  cylindrical  shells  under  impulsive  radial  pressure  [4.2]  and  in  thin 
bars  under  axial  impact  [4.  3].  In  the.e  latter  problem.,  the  \c i>  »hc*t 
wavelengths  can  be  attributed  directly  to  the  very  high  thrusts  compared 
to  those  in  static  buckling.  In  cylindrical  shells  under  dynamic  axial 
thrust,  however,  the  thrust  does  not  have  to  be  greater  than  the  classical 
static  buckling  load  to  produce  wavelengths  shorter  than  those  observed 
in  large  deflection  static  buckling.  Using  high-speed  photography  of  a 
shell  buckling  elastically  under  static  thrust,  Almroth,  Holmes,  and 
Brush  [4.4]  showed  that  the  shell  initially  deformed  into  diamonds  of 
roughly  half  the  size  of  those  in  the  final  buckled  shape.  In  other  experi¬ 
ments,  by  placing  a  rigid  mandrel  inside  the  shell  to  limit  the  amplitude 
of  post-buckling  deformation,  they  were  able  to  keep  the  shell  buckled 
in  the  shorter  wavelength  pattern,  Thus,  in  a  cylindrical  shell  under 
dynamic  axial  thrust,  the  modes  of  deformation  can  be  expected  to  depend 
in  general  upon  the  amplitude  of  deformation  as  well  as  upon  the  magni¬ 
tude  of  the  thrust. 


Th-se  complications  in  th«  cylindrical  shell  have  led  theoretical 
investigations  of  dynamic  buckling  toward  a  large  deflection  theory  as 
was  found  to  be  required  for  static  buckling.  The  general  procedure 
has  been  to  assume  axial  and  circumferential  buckle  wavelengths  and  to 
numerically  integrate  equations  derived  from  large  deflection  shell 
theory.  The  coefficients  of  a  few  terms  giving  roughly  the  diamond 
buckle  shape  are  taken  as  generalized  coordinates  and  the  Galerkin  pro¬ 
cedure  is  used  to  derive  equations  governing  these  coefficients.  Buckling 
is  said  to  occur  when  the  deflections  increase  abruptly  with  load  or  end 
shortening.  The  procedure  is  repeated  over  a  range  of  axial  and  cir¬ 
cumferential  wavelengths  to  find  the  combination  that  gives  the  lowest 
buckling  load.  Coppa  and  Nash  [4.  5]used  a  twodegree -cf -freedom  system, 
taking  the  end  shortening  to  increase  linearly  with  time.  Roth  and 
I»lcsr.cr[4  use  J  a  four  degree -of -freedom  system  and  a  step  pressure. 

No  comparison  of  theae  theories  with  experiment  was  made. 

The  drawback  of  the  above  procedure  is  that  two  of  the  essential 
unknowns  of  the  problem,  the  axial  and  circumferential  wavelengths  of 
the  buckles,  must  be  assumed  at  the  outset  of  the  integration.  To 
determine  the  wavelengths  into  which  the  shell  actually  buckles,  the 
procedure  must  be  repeated  over  and  over  for  many  combinations  of 
axial  and  circumferential  wavelengths.  For  buckling  in  which  the  impact 
stress  is  greater  than  the  static  classical  buckling  stress,  this  short¬ 
coming  is  overcome  in  the  present  paper  by  using  linear  theory  to  study 
the  wave  selection  process  in  the  early  motion  while  the  Jeflections  are 
still  small.  This  allows  application  of  the  analytical  procedure  used  in 
the  beam  under  axial  impact  [4.  3], The  buckled  form  is  not  specified 
a  priori;  instead  the  motion  is  expanded  into  the  linear  vibration  modes 
of  the  shell  so  that  it  is  free  to  deform  into  its  own  "preferred"  shape. 
Initial  imperfections  are  taken  in  the  form  of  white  noise  and  a  statistical 
analyses  gives  the  expected  values  for  the  axial  and  circumferential 
wavelengths.  If  the  initial  imperfections  are  sufficiently  small,  the  ratio 
of  amplitudes  of  the  preferred  modes  to  the  amplitudes  of  the  imperfec¬ 
tions  can  become  quite  large  within  the  applicability  of  small  deflection 
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theory,  and  this  dominance  of  the  preferred  linear  wavelengths  could 
carry  over  into  the  large  deflection  motion. 

Experiments  are  presented  which  give  final  buckled  forms  very 
close  to  those  predicted  by  the  theory.  Very  high-speed  motion  pictures 
(240,  000  frames  per  sec)  show  buckles  forming  during  the  elastic  and 
plastic  motion  and  demonstrate  that  the  time  scale  of  the  buckle  for¬ 
mation  also  agrees  with  the  theory.  The  axial  wavelengths  are  much 
shorter  than  in  static  buckling  and  the  observed  aspect  ratios  (circum¬ 
ferential  to  axial)  are  much  larger,  averaging  about  3  compared  to  1  in 
static  buckling.  This  varied  considerably  and  in  many  experiments 
nearly  symmetric  (very  large  aspect  ratio)  buckling  was  observed. 

4.  2  RESPONSE  Ob'  THE  NORMAL  MODES 

In  the  experiments,  cylindrical  shells  are  impacted  at  one  end  by 
a  massive  ring  which  sends  a  step  axial  stress  wave  down  the  shell. 

T'ne  shell  buckles  before  any  signal  is  received  from  the  opposite  (free) 
end  so  that  in  the  theory  the  shell  is  taken  arbitrarily  long.  Also,  the 
thrust  is  assumed  uniform  throughout  the  length  of  the  shell,  thus  neglect* 
ing  any  effects  of  the  moving  axial  stress  front.  This  was  shown  to  be  a 
reasonable  assumption  for  the  axially  impacted  strip  [4.  3],  and  appears 
to  be  valid  here  also,  since  as  shown  below,  the  axial  wavelengths  and 
formation  times  of  the  buckles  in  the  strip  and  cylinder  are  quite  similar. 
Thus,  we  consider  a  long,  thin  cylindrical  shell  subjected  to  a  suddenly 
applied  constant  compressive  force. 

To  study  the  small  early  displacements,  it  is  sufficient  to  use  a 
small  deflection  linear  theory.  It  can  be  further  shown  that  simplified 
Donnell-type  equations  give  satisfactory  accuracy  for  the  present  problem. 


If  more  complete  equations,  such  as  given  in  [4.  8],  are  used, the  final 
equation  (12)  here  results  if  terms  multiplying  the  axial  strain  c  *  N/Eh 
and  terms  containing  (h/a)^  are  omitted  as  small  compared  to  terms 
without  these  small  multipliers. 


53 


Denoting  time  by  t  and  axial  force  per  unit  arc  length  by  N  ,  and 
taking,  as  shown  in  Fig.  4.  1,  x  and  y  to  be  the  axial  and  circumferential 


m  •  NUMBER  OF  AXIAL  HALF-WAVES 

n  •  number  of  circumferential  waves 


FIG.  4.1  COORDINATES  AND  SHELL  NOMENCLATURE 

coordinates,  and  w  to  be  the  radially  inward  displacement  measured 
from  the  initial  perturbed  displacement  wq  ,  these  equations  become 

d2  ,  ,  .  ,  .  h  d2F  , . 

DV  ^  i  N  1  n  (w  +  W  )  +  pn  i,,i,,,|y  •  —  - mn  c  0  (4*1) 

*x£  °  ar  4  a*2 

V4F  *  7  ^7  s  0  <4*2> 

*  a*£ 

where  F  is  Airy's  stress  function  for  mid-surface  stresses  produced 

by  the  buckling,  a  is  the  shell  radius,  h  is  the  thickness,  c  is  the 

density  ,  D  =  Eh^/12(l-v^)  is  the  flexural  stiffness  and  V*  is  the 
2  2  2  2  2 

operator  (d  /dx  ♦  a  /dy  1  .  These  are  equations  (3)  and  (4)  from 
(4.7)  with  terms  appropriate  to  the  present  problem  retained. 
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Using  non-dimensional  variables 


xJb  ' 


N 


V  zhD~ 


[‘ i ■  3) 


equations  (4.  1)  and  (4.2)  become 
4  3^ 

"  W  +  - y  (w  +  W  )  r  W 


h  32F 

^  7^ 


=  o 


(4.4) 


r-4  ED  S  w  ~ 

v  F  +  — r-r  =  0 


aN 


a?' 


(4.5) 


where  the  operator  is  now  (3^/35^  +  c^/3r^)^  .  and  (*)  =  3/3* 


To  simplify  the  mathematics,  the  shell  will  be  considered  simply 
supported.  The  edge  condition  at  the  impacted  end  of  the  shells  in  the 
experiments  is  more  closely  approximated  by  a  clamped  boundary,  but 
the  experiments  indicate  that  this  boundary  condition  does  nut  signifi¬ 
cantly  affect  the  buckling.  The  conditions  of  simple  support  are  satisfied 
by 


w 


(?. 


00 


T)  = 


■El 

m  ~  1  n  - 1 


W  {*)  sin  q  *  sin  8  r 
mn  m  n 


(4.  6} 


F  {?  . 


t) 


m=d 


F  (-)  sin  a  faint?  *• 

mn  m  -  n 


(4.  7) 


who  re 


n 


n  fp" 

a  4N 


(4.  H) 


and  L  is  the  length  of  the  shell.  The  initial  displacement  w  is  also 
expanded  into  a  Fourier  series: 


w  (' 


■tt 

m  =  1  n  =  1 


a  sin 
mn 


3  9  sin 

n> ' 


(4.9) 
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Equation  (4.  7)  is  now  substituted  into  (4.5)  and  the  result  is  substituted 
with  (4.  6)  into  (4.  4)  to  give  the  following  equation  for  the  amplitudes 


W  of  the  normal  modes, 
mn 


,  r,-2  .  „2*2  2  ,  EDh  am  t  _  2 

W  +  >•  ^am  +  - 7  — 1 - — 7*" 7  J 

mn  m  n  m  £.,<•  ,  £  ,  nc  .&  mn  m  mn 

a  N  (a  +8) 
m  n 


(4.  10) 


In  the  limit  as  L  — »*oo  i  a  becomes  a  continuous  variable  and,  for 

m 

convenience,  3  will  also  be  treated  as  a  continuous  variable  so  that 
n 

in  the  following  Wmn<T  )  is  replaced  by  W(a,  0,  t)  and  the  subscripts 
on  a  and  0  are  dropped.  Noting  further  that 


EDh 

~T~J 
a  N 


^>2 


where 


ct 


V  3(l-v2) 


h 

a 


(4.11) 


is  the  classical  static  buckling  stress,  (4.  10)  can  be  written 


g  +  k(a,  0  )g  =  a 


(4. 12) 


where 


k(a,  8)  =  (a2+02)2  -  a2  +  t  (  ~-)2  -^~TT 

4  a  (a2+  02)2 


(4. 13) 


and 


_/.  «  _  *  _  W(a,  8  »  t  ) 

g(cx,  8,  t)  -  -Ya-;--py- 


(4.  14) 


is  the  amplification  function  or  "filter"  characteristic  of  the  shell  under 
a  given  axial  stress  a  =  N/h  . 
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;ind  the  hyperbolic  form  is  taken  for  k(a,  6  )  <  0.  Figure  4.2  gives 
example  pkrs  of  g(a  ,  0,  t)  for  t  =■  4  and  t  =  8  .  The  dependence 
of  these  amplitudes  on  the  axial  wave  number  a  is  similar  tc  that  in 
the  buckling  bar,  exhibiting  a  pronounced  hump  of  "preferred"  modes 
near  a  =  1/^2  .  The  circumferential  wave  number  of  the  most  ampli¬ 
fied  mode  is  0  =  0  ,  i. e. ,  a  symmetric  mode,  but  there  is  an  appreciable 
bandwidth  of  amplified  modes  in  both  the  axial  and  circumferential 
directions.  Comparing  the  curves  for  r  =  4  and  t  =  8  we  see  that  as 
the  motion  proceeds,  the  bandwidth  in  both  directions  decreases,  tending 
in  the  limit  to  produce  a  fixed  axial  wavelength  corresponding  to  the 
wavenumber  a  «  l/'Jz  ,  but  tending  toward  longer  wavelengths  in  the 
circumferential  direction,  ultimately  approaching  a  simple  symmetric 
pattern  corresponding  to  0  =  0  (i.e. ,  n  =  0).  If  the  imperfections  are 
large  enough,  however,  non-linear  effects  would  dominate  before  the 
symmetric  mode  is  attained,  producing  final  buckles  with  a  finite  wave¬ 
length  in  both  directions.  Buckling  of  both  types  were  observed  in  the 
experiments  described  below. 


From  the  definition  of  k(a(  0  )  in  (4.  13)  we  see  that  the  normalized 

amplification  curves  g(a,  0,  t)  depend  only  on  the  ratio  o/o^  of  the 

applied  axial  stress  to  the  classical  static  buckling  stress.  Larger 

values  give  faster  growth  and  narrower  bandwidths,  but  the  general 

shape  of  the  amplification  curves  remains  about  the  same  for  any  o/oc^>  1. 

In  the  next  section  it  will  be  shown  that  the  expected  values  of  the  buckle 

rn —  f— 

wavelengths  are  very  nearly  proportional  to ^  =  r/«y  e  ,  where 
r  =  h/<\/l2(l-v^)  and  e  =  o/E  is  the  strain  from  the  axial  thrust. 
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0C49M-2 


FIG.  4'. 2  AMPLIFICATION  FUNCTION  g  vs.  AXIAI.  AND  CIRCUMFERENTIAL 
WAVE  NUMBERS  a  AND  fi  FOR  Jocl  -  2.5 


Figure  4.  3  gives  a  plot  of  the  maximum  value  of  g  (at  the  peak  of 
the  hump  in  Fig.  4.2)  plotted  against  time  for  several  values  of  o/o  . 

C  v 

As  a/h  — ■*■»  we  see  from  (4.  11)  that  a  .  — *»  0  so  that  the  curve  for 

Cv 


c/a  — ►  CO  is  that  for  a  buckling  plate  (or  rectangular  bar  if  the  factor 

C  v 

l-v^  is  omitted).  For  a/a^^>  2  the  time  required  to  attain  large 
amplification  is  only  slightly  greater  than  for  the  plate. 


FIG.  4.3  GROWTH  OF  MAXIMUM  AMPLIFICATION  WITH  TIME 

4.  3  BUCKLING  FROM  RANDOM  IMPERFECTIONS 

To  compare  theory  with  experiment,  it  is  necessary  to  assume  some 
form  of  imperfections  in  order  to  specify  a(a,  0  )  and  hence  to  compute 
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W(a,  8*  t  )  from  (4.  14).  Experiments  on  impact  buckling  of  bars  in 

[4.  3]  demonstrated  that  a  good  description  of  the  observed  buckles  was 

obtained  by  assuming  imperfections  in  the  form  of  white  noise.  These 

allow  the  buckle  wavelengths  to  be  dictated  by  the  magnitude  of  the  thrust, 

giving  a  random  scatter  of  wavelengths  with  a  mean  and  standard  deviation 

both  inversely  proportional  to  the  square  root  of  the  thrust,  in  agreement 

with  experimental  observations.  It  is  expected  that  such  an  assumption 

will  also  be  reasonable  in  the  present  problem.  Thus,  we  assume  the 

a(a,  0  )  are  random  normal  with  zero  mean  and  constant  variance  over 

all  a  and  0  in  the  amplified  band  of  interest.  With  this  assumption, 

the  power  spectral  density  of  the  modal  amplitudes  W(a,  8  ,  T  )  is 
2 

proportional  to  g  (a,  0  ,  t  ). 


Having  the  power  spectral  density,  the  statistics  of  the  buckled 
shape  can  be  computed.  For  the  bar  [4.  3],  complete  statistics  of  the 
wavelengths  were  computed  using  a  Monte  Carlo  technique.  The  simpler 
problem  of  determining  the  mean  wavelengths  can  be  solved  analytically 
and  will  be  undertaken  here  for  the  shell.  From  Rice  [4.  9]  the  expected 
number  of  zeros  in  an  interval  (Sj  ,  s^)  of  a  random  function  f  *  F(s) 
is 

*2 


Op 

j  <*■  J  |q|  p(0.  q;  »X*q 
a.  -oo 


(4.  16) 


where  P(f  ,  q  ;  s)  is  the  probability  density  function  for  the  variables 


i  •  FIs) 
_  .  dF 

11  ■  a? 
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V 


For  our  problem  we  have  then,  that  the  expected  number  of  zeros  in 
the  §  direction  in  the  interval  (0  ,  l  ),  for  a  specified  value  of  n  and 


t  00 

R  =  J  d§  f  Z  P(0.  z;?,  n,  T)dz 
0-00  * 


(4.  17) 


/n* 

where  t  =  Lyg  and  P{w  ,  z;  ?,  ri  ,  t)  is  the  probability  density 
function  of  the  variables 

W  =  W(?  ,  T]  ,  t  ) 


and  ?  ,  ri  ,  t  are  carried  along  as  parameters.  Since  it  has  been 
assumed  that  the  initial  imperfections  have  a  Gaussian  probability 
distribution  (or,  alternately,  if  we  apply  the  central  limit  theorem  for 
more  general  imperfection  statistics)  then  the  linal  buckled  form  will 
also  have  a  Gaussian  distribution.  Thus,  the  distribution  function  has 
the  form 


P(w,  *;  ?  ,  r  ,  t  )  = 


2n  ^11^22  *  ^12  > 


where 


~*n*2  *  2^i2W8  m*zz*Z 
2  ^  11^22  “  ^12^ 


»  <w‘> 

Hz  s  <,2> 

P12  »  <WS> 


and  <>  indicates  ensemble  average. 


(4.18) 


-  81 


Our  results  are  simplified  if  it  is  assumed  that  the  imperfections 
and  final  buckled  form  are  stationary  (in  space)  sc  that  in  place  of  (4.  6) 

we  Cum  write 


w(S,  r\.  t)  ■ 


j  £  Wm»<T>  ■*"'«««  '  ®m>  *in  <V>-  «„>  «•  19> 

m=  I  n*  1 

where  the  and  8  are  each  uniformly  distributed  over  the  interval 
m  n 

(0,  2tt).  With  this  assumption  the  statistics  of  the  buckled  form  are 
independent  of  5  and  ri  ,  which  we  would  expect  to  be  justifiable  for 
waves  at  a  sufficient  distance  from  the  ends  of  tie  cylinder.  In  fact, 
the  results  of  a  Monte  Carlo  computation  including  end  effects  for  the 
bar  [4.  3],  which  is  similar  in  form  to  the  cylinder,  indicate  satisfactory 

agreement  with  the  stationary  process  assumption  even  for  the  first 
wave  from  the  simple  support. 


Using  (4. 19)  it  can  be  easily  shown  that 

co  oo  oo  ao 

u  a  <w2>  *  JJ  W2(a.  8.  *r)dad0  =  c  2f j  g2(a.  8.  t  )dad8 
0  0  0  0  (4. 


20a) 


U22  *  <«  > 


J  £  a2w2(a,  8,  t)dad8  *ol  J j  a2g2(a. 


8,  T)dadB 

(4.20b) 


*  <»*  >  ■  0 


(4.20c) 


where 


2  2 

9  *  <  a  (a ,  8  )  >  *  constant 


and  the  sums  have  been  replaced  by  integrals. 


Making  use  of  =  ®  from  (4.20c)  in  (4.  18)  and  substituting  the 
resulting  expression  in  (4.  17)  gives  the  expected  number  of  zeros  in  the 
length  l  , 


R  =  —  (  — )  1/2 


(4.21) 


The  mean  wavelength,  as  measured  between  alternate  zero  crossings 
is  simply 


=  W  =  2n(r1)1/2 
5  «  ^22 


(4.22) 


or,  upon  using  (4.20) 


=  *‘T 


oo  ac 

l 


g  (a ,  8 ,  t  )dadg 


J  J  a2g2(a,  S,  T)dad0 


(4.23) 


The  mean  wavelength  in  the  circumferential  direction  can  similarly  be 
shown  to  be 


X  =  2tt 
T1 


op  00  2 

J  J  g  (a ,  8  ,  T  )dad8 
0  0 _ 

}h  g2(a  .  8  ,  t  )dadB 


(4.24) 


Equations  (4.23)  and  (4.24)  were  integrated  numerically  over  the 
area  0&s&2,0*8*2  of  significant  amplification  and  the  results 
are  presented  in  Fig.  4.4.  In  Fig.  4.4(a)  the  variation  of  X^  and  X^ 
with  t  is  shown.  Bands  have  been  drawn  which  contain  the  curve*  for 

o to  .  in  the  range  from  1.  12  to  4.  0.  The  narrowness  of  these  bands 

ct 

indicates  that  the  most  significant  parameters  affecting  wavelength  are 
those  given  in  the  normalisation  in  (4.  3). 
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ASPECT  RATIO,  X,/  X<  NORMALIZED  WAVELENGTHS.  AMO 


FIG.  4.4  AXIAL  AND  CIRCUMFERENTIAL  MEAN  WAVELENGTHS  v».  TIME 
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A'so  fro;n  Fig.  4.  4(a),  we  see  that  the  mean  axial  wavelength 
increases  only  aiightly  with  t  for  t  >  4  and  quickly  approaches  a 
"preferred"  wavelength.  For  ratios  of  c/c^  nearing  unity  the  preferred 
wavelength  becomes  the  classical  static  wavelength  >.p  =  Z»Jz  The  mean 
circumferential  wavelength,  however,  increases  with  t  without  approaching 
an  asymtotic  value.  Thus,  as  discussed  in  the  previous  section,  the  mean 
cixc.un.fv.iCr*"1  "•■"’“lengths  actual'v  observed  it  largo  deflection  buckling 
will  probably  depend  on  the  magnitude  of  the  initial  imperfections,  which 
determines  the  duration  for  which  this  small  deflection  theory  is  applicable. 
Smaller  i  nperfections  can  grow  for  a  longer  time  and  from  Fig.  4.4(a) 
we  would  expect  to  see  longer  circumferential  wavelengths. 

Figure  4.4(b)  gives  a  plot  of  aspect  ratio  !\f  .  If  nonlinear 
effects  begin  to  dominate  at,  say  -  =  7  ,  and  subsequent  buckling  proceeds 

with  a  fixed  pattern  (one  would  expect  a  fixed  pattern  to  be  eventually 
established,  as  confirmed  in  the  experiments),  Fig.  4.4(b)  indicates 
that  the  aspect  ratio  would  be  about  3.  3.  This  is  discussed  again  in 
reference  to  the  experiments. 

4.4  EXPERIMENTS 

A  small-deflection  linear  theory  applied  to  the  static  buckling  of 
cylindrical  shells  under  end  load  is  notoriously  inadequate  to  predict 
experimental  buckling  loads  except  for  very  nearly  perfect  shells.  Also, 
the. shells  ultimately  buckle  into  a  form  very  different  from  the  buckling 
mode  of  the  small- deflection  theory.  It  is  not  obvious,  then,  that  the 
simple  linear  theory  given  here  should  reasonably  predict  the  large 
deflection  dynamic  buckling  of  such  a  shell.  It  was  thought  however, 
that  the  small-deflection  theory  would  give  promise  for  the  dynamic 
problem  because,  once  started  in  the  linear  theory  pattern,  the  shell 
would  continue  to  deform  in  this  pattern  and  not  have  time  to  convert  to 
another  pattern.  The  experiments  described  here  were  run  simultaneously 
with  the  theoretical  investigation  to  determine  whether  such  a  theory  should 
be  pursued  at  all. 
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Previous  expe r in iento rs  reported  only  the  final  shapes  (diamond 
buckles)  after  very  severe  and  prolonged  buckling,  and  showed  high¬ 
speed  (about  15,  000  frames/secl  motion  pictures  whicn  were  neverthelr -a 
at  a  speed  too  slow  to  show  the  early  buckling  process.  In  order  to  record 
a  y  motion,  tne  present  experimental  arrangement,  shown  in  Fig.  4.  b, 
incorporated  three  u..-iue  features:  (1)  the  shells  were  free  at  the  end 
opposite  the  impact  so  that  the  compressive  impact  stress  would  have  a 
duration  (at  the  impacted  end)  of, at  most, equal  to  the  transit  time  2L/c 
of  the  longitudinal  stress  wave  up  and  down  the  shel1  This  allowed  termi¬ 
nal  observation  of  early  permanent  buckling  not  obliterated  by  later  folding. 
(2)  The  impacted  end  was  bonded  with  epoxy  cement  to  a  relatively  massive 

inside  ring,  and  to  a  thin  outside  ring  to  provide  a  clamped  boundary  to 
the  shell.  Without  these  rings  severe  crimping  at  the  end  rapidly  lowered 
the  thrust  and  made  comparison  to  theory  impossible.  (3)  The  massive 
end  ring  was  accelerated  explosively  so  that  the  time  and  simultaneity  of 
impact  could  be  controlled  to  within  about  2  microseconds.  This  allowed 
the  use  of  a  Beckman- Whitley  framing  camera  running  at  240,  000 
frames /sec,  fast  enough  to  see  the  details  of  early  wave  formation. 

Test  shells  were  made  from  0.  002 7-inch-thick  5052-H19 
aluminum  «  heet  rolled  to  a  three -inch  diameter  with  a  lapped  seam 
held  with  cloth  tape.  The  steel  ring  at  the  lower  end  of  each  shell 
served  as  the  ” impacting”  mass  and  weighed  12  times  the  weight  of  the 
shell  so  that  its  change  in  velocity  during  the  impact,  and  hence  the 
change  in  impact  stress,  was  small.  The  ring  was  accelerated  by  placing 
it  on  the  heavy  steel  anvil  bar  and  detonating  a  sheet  explosive  charge  on 
the  opposite  end  of  the  anvil.  The  explosive  sent  a  steep-fronted  shock 
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TEST  CYLINDER,  EROM  5052 -MI9 
aluminum  sheet 


Fig  4  5  EXPERIMENTAL  ARRANGEMENT 


into  the  anvil  which  entered  the  ring  and  bcunced  it.  off  in  much  the  same 
Vv ay  ?s  are  the  end  pellets  from  a  Hopkinson  bar  [  4.  10} .  The  pressure 
gradient  behind  the  shock  was  sufficiently  shallow  that  the  reverberating 
stresses  in  the  ring  were  small  and  it  is  assumed  that  the  ring  was 
stress -free  when  it  "impacted"  the  cylinder.  To  insure  good  transmission 
of  the  shock  from  the  anvil  to  the  ring,  the  contact  surfaces  were  lapped. 

The  ring  velocity  from  various  thicknesses  of  sheet  explosive  was 
determined  in  a  separate  series  of  experiments.  Since  the  ring  is  very 
massive  compared  to  the  shell,  the  impact  stress  in  the  shell  is  assumed 
to  be  pcv  ,  where  pc  is  the  acoustic  impedance  of  aluminum  shell,  and 
v  is  the  velocity  of  the  steel  ring.  The  axial  strain  corresponding  to 
this  stress  is 


where  c  is  the  axial  wave  speed  in  the  shell. 

Figure  4.  6  shows  a  sequence  of  framing  camera  pictures  taken  in 
an  experiment  in  which  v  =  340  in/sec  which,  with  c  =  200,  000  in/sec, 
gives  e  =  0.0017.  Substituting  a/h  =  1.50/0.0027  =  555  into  (4.11) 
yields  e  =  0.00109  so  that  e/e  -  1.56.  Thus,  the  impact  stress 

C  v  C-C 

is  1.56  times  the  classical  static  buckling  3tres9.  Only  alternate  frames 
from  the  original  record  are  shown,  giving  8.  33  psec  between  frames 
and  an  exposure  time  of  about  1.4  p sec.  Normalized  time  t  can  be 
computed  from  (4.  3)  which  yields  t  =  (ce/r)t  =  0.  38t  .  On  the  original 
record,  very  tiny  a.splacerr.ents  could  be  discerned  at  t  =  7  psec  (  t  =  2.  6). 
At  t  =  1 1 . 2  psec  (t  =  4.  2)  in  Fig.  4.6,  small  wrinkles  near  the  bottom 
of  the  shell  are  clearly  visible.  By  t  =  27.  8  psec  *t  =  10.  6)  these  wrinkles 
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FIG.  4.6  AXIALLY  IMPACTED  CYLINDRICAL  SHELL 
Time  it  from  initial  impact  at  lower  end 


are  definitely  taking  on  the  diamond  pattern,  indicating  that  nonlinear 
effects  are  predominating.  Thus,  from  t  =  2.  6  to  t=  10.6  the  displace¬ 
ments  grow  from  being  just  visible  to  amplitudes  so  large  that  nonlinear 
effects  dominate.  This  agrees  very  well  with  the  period  of  first  rapid 
growth  given  by  the  theory.  From  Fig.  4.3  at  ala  =  1.  5  we  see  that 

the  amplification  grows  from  g  =  3  at  t=  2.6  to  g  =  100  at 

mix  max 

T  =  10.  6  .  It  is  probable,  as  discussed  later,  that  nonlinear  effects  begin 
to  dominate  at  an  intermediate  time  of  about  t  =  7  . 

At  frame  t  =  27.  8  psec  ,  the  beginning  of  buckling  farther  up  the 
shell  is  evident.  In  the  following  frames  these  buckles  stay  fixed  in 
position  and  grow  in  amplitude.  Other  buckles  appear  at  areas  in  between, 
at  which  initial  imperfections  were  probably  smaller.  The  tiny  ripples 
just  above  the  buckles  at  the  lower  end  characteristically  appeared  in  all 
of  the  experiments.  These  are  most  visible  at  t  =  44.  5  psec. 

To  compare  the  wavelengths  in  Fig.  4.  6  with  theory,  the  expected 
value  of  the  axial  wavelength  is  determined  by  substituting  X^  =  8.  9 
from  Fig.  4.4  into  (4.3),  which  yields 


Using  c  =  0.  0017  and  h  =  0.  0027  inch  gives  X^  =  0.  18  inch.  From 
Fig.  4.  6,  the  average  length  from  15  waves  is  X^  =  0.20  inch,  only 
10%  greater  than  the  theoretical  mean. 

The  permanent  buckles  remaining  in  several  shells  after  impact 
are  shown  in  Fig.  4.  7.  Figure  4.  7(a)  is  of  the  same  shell  as  in  Fig. 

4.  6  and  it  is  apparent  that  all  but  the  lower  few  buckles  in  Fig.  4.  6  were 
elastic.  Most  of  the  kinks  away  from  the  end  in  Fig.  4.  7(a)  were  the 
result  of  difficulties  in  retrieving  the  shell  (it  was  projected  12  feet  into 
the  air).  In  general,  permanent  buckles  were  confined  to  the  area  close 
to  the  impacted  end  as  seen  in  all  the  shells  in  Fig.  4.  7.  The  most 
striking  feature  of  these  buckles  is  their  very  small  size  and  large  aspect 
ratio  as  compared  to  static  buckles.  For  comparison,  Fig.  4.8  shows  a 
statically  buckled  shell  of  the  same  material,  wall  thickness  and  diameter 
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FIG.  4.7  PERMANENT  BUCKLES  FROM 
AXIAL  IMPACT  (oppotito  and 
wot  frt«) 
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FIG.  4.8  STATIC  BUCKLING  OF 
N  IDENTICAL  SHELL 


as  those  in  Figs.  4.6  and  4.7.  The  circumferential  and  axial  wave¬ 
lengths  are  about  1.  3  inches,  giving  an  aspect  of  unity.  The  dynamic 
buckles  in  Fig.  4.  7  have  a  mean  aspect  ratio  of  3.  3  and  axial  wave¬ 
lengths  of  only  0.  18  inch.  Several  shells  are  shown  to  illustrate  the 
wide  range  of  aspect  ratio  observed.  There  is  little  doubt  that  the 
shape  of  the  buckles  was  strongly  influenced  by  large  elastic  and  plastic 
deformations,  but  it  is  significant  that  the  large  aspect  ratio  is  suggested 
by  the  linear  theory  of  the  preceding  sections. 

Circumferential  wavelengths  were  measured  from  the  permanent 
buckles  in  these  and  two  other  shells  and  the  results  are  summarized 
in  Fig.  4.9.  Aspect  ratio  is  plotted  rather  than  wavelength  to  empha¬ 
size  the  difference  between  these  buckled  forms  and  static  buckles.  The 
values  range  from  2.2  to  7.  1  with  a  mean  of  3.3.  These  are  much  larger 
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FIG.  4.9  HISTOGRAM  OF  OBSERVED  ASPECT  RATIOS 


than  the  static  ratio  of  about  1.  It  is  difficult  to  compare  these  wave¬ 
lengths  with  the  linear  theory  because,  as  shown  in  Fig.  4.4,  the  theo¬ 
retical  wavelengths  continually  increase  with  time.  However,  from 
Fig.  4.4(b),  the  ratio  of  the  theoretical  average  circumferential  wave¬ 
length  to  the  average  axial  wavelength  is  3.  3  at  t  =  7  .  Also,  the  photo¬ 
graphs  in  Fig.  4.  6  show  that  this  is  about  the  time  at  which  the  buckling 
amplitude  becomes  large. 

This  excellent  agreement  between  theoretical  and  observed  axial 
wavelengths,  aspect  ratios,  and  buckling  times,  confirms  our  original 
premise  that  early  small  deflection  buckling  dictates  the  pattern  into 
which  large  deflection  buckles  form.  If  more  details  of  the  buckling  are 
desired,  displacements  in  the  "preferred"  wavelengths  from  the  linear 
theory  could  be  used  as  initial  conditions  in  a  large  deflection  theory. 

For  calculating  the  threshold  of  buckling  from  axial  impact,  it  may  be 
sufficient  to  merely  define  some  magnitude  of  the  linear  amplification 
function  as  indicative  of  buckling. 

Crumpling  of  cylindrical  shells  used  as  an  energy  absorbing 
mechanism  presents  a  much  more  complicated  theoretical  problem  in 
which  the  final  b  ickling  may  be  in  a  mode  much  different  from  that 
observed  here.  This  occurs  because  deformation  in  the  very  short 
axial  wavelengths  discussed  here  becomes  so  large  that  severe  crumpling 


occurs  at  one  or  both  ends  (the  duration  of  the  initial  thrust  is  longest 
at  one  end  because  of  the  finite  speed  of  axial  propagation).  This 
greatly  reduces  the  thrust  and  subsequent  buckling  in  the  remainder  of 
the  shell  probably  occurs  at  a  fluctuating  thrust  which  depends  on  the 
buckling  itself.  An  example  of  such  buckling  is  given  in  Fig.  4.  10. 


FIG.  4.10  BUCKLING  WITH  A  HEAVY  MASS 
ON  BOTH  ENDS,  a/acl  =  1.56 
AT  INITIAL  IMPACT 


The  impact  velocity  and  shell  parameters  are  the  same  as  in  Fig.  4.  6, 
but  now  rigid  masses  are  attached  at  both  ends  so  that  buckling  pro¬ 
ceeds  for  a  longer  time.  Buckling  away  from  the  ends  bears  little 
resemblance  to  the  impact  buckling  in  Figs.  4.  6  and  4. 7,  and  in  fact, 
is  strikingly  similar  to  the  static  buckling  in  Fig.  4.  8.  However,  the 
crimping  near  each  end  in  Fig.  4.  10  is  of  the  same  form  as  the  impact 
buckling  in  Figs.  4.  6  and  4.  7.  Thus,  the  present  theory  predicts  the 
early  phase  of  buckling  and  suggests  that  away  from  the  ends  the 
transition  to  the  final  buckled  state  is  very  complex  and  can  be  studied 
only  by  including  the  interaction  between  the  buckling  and  axial  thrust. 
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5.  DYNAMIC  BUCKLING  OF  A  CYLINDRICAL  SHELL 
OF  A  STRAIN-RATE  SENSITIVE  MATERIAL 


5.  1  INTRODUCTION 

In  this  section  the  dynamic  flow  buckling  of  a  cylinder  subjected 
to  a  uniform  inward  radial  impulse  for  the  case  of  a  material  which 
does  not  exhibit  strain  hardening  but  is  strain-rate  sensitive  is  investi¬ 
gated.  Cylinders  of  fully  annealed  1015  steel,  a  material  known  to  be 
strain-rate  sensitive  [5.1,  5.2,  5.3,  5.4],  were  subjected  to  impul¬ 
sive  loads  which  produced  final  wrinkled  forms  similar  to  that  shown 
in  Fig.  5.  1  with  average  wavelengths  that  were  quite  reproducible. 


FIG.  5.1  TYPICAL  BUCKLED  SHAPE  OF  CIRCULAR  CYLINDRICAL  SHELL 
OF  FULLY  ANNEALED  1015  STEEL  —  CYLINDER  NO.  lb 
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These  wavelengths  are  compared  with  the  predictions  of  a  theory  of 
flow  buckling  which  assumes  that  the  material  is  linearly  visco-plastic. 
The  reasonably  good  correlation  suggests  that  this  behavior  could 
account  for  the  formation  of  the  final  buckled  shape. 

Abrahamson  and  Goodier  [  5.  5]  proposed  a  theory  of  the  for¬ 
mation  of  the  wrinkled  shape  based  on  the  dynamics  of  a  metal  shell 
which  deforms  plastically  with  linear  strain  hardening.  This  latter 
property  is  crucial  to  the  theory  because  without  it  divergent  series 
result.  Experimental  and  predicted  numbers  of  waves  were  in  good 
agreement.  Goodier  and  Mclvor  [  5.  6]  gave  a  similar  theoretical 
treatment  for  a  material  which  exhibited  both  linear  strain  hardening 
and  linear  strain-rate  sensitivity.  In  particular,  theoretical  results 
are  given  for  linear  strain-rate  dependency  with  no  strain  hardening 
(visco-plastic),  which  is  the  behavior  assumed  here.  An  outline  of 
the  derivation  of  the  equations  of  motion  and  their  solutions  is  given 
later. 


5.2  DYNAMIC  FLOW  BUCKLING 

When  the  elements  of  a  perfect  cylindrical  shell  simultaneously 
receive  the  same  large  initial  inward  radial  velocity  the  material  flows 
into  a  uniform  cylindrical  shell  of  smaller  radius  and  thicker  wall  until 
the  initial  kinetic  energy  has  been  converted  into  plastic  work.  Inevitable 
small  imperfections,  such  as  in  the  uniformity  of  the  initial  velocity, 
cause  perturbations  from  this  uniform  converging  motion.  Departures 
from  the  circular  form  are  amplified  by  the  action  of  the  compressive 
circumferential  stress.  These  departures  are  considered  small  enough 
to  regard  the  compressive  circumferential  strain  as  increasing  through¬ 
out  the  motion.  Thus  at  two  points  A  and  B,  on  the  outside  and  inside 
respectively,  as  shown  in  Fig.  5.2,  the  circumferential  strain  increases 
but  more  rapidly  at  B  than  A  due  to  the  slight  curvature  perturbation. 
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Or>  the  stress  vs.  strain-rate 
diagram  in  Fig.  5.2  the  representa¬ 
tive  points  A  '  and  B  are  both 
moving  to  the  left  with  B  ‘  lagging 
behind  A  ; .  The  strain-rate  law  is 
taken  as 

c  =  o  *  pc,  (5.  1) 

o  t  ' 

where  c  and  et  are  the  stress 
and  the  strain  rate  (subscripts  t 
denote  partial  differentiation  with 
respect  to  time)  and  p  i3  the  slope 
of  the  straight  line.  Owing  to  the  difference  of  circumferential  stress 
between  A*  and  B  ‘  there  is  a  bending  moment  M  in  the  section  AB 
of  the  shell.  For  a  unit  axial  length  of  cylinder  regarded  as  a  ring, 
integrating  (5.  1)  over  the  cross  section,  with  the  assumption  that  plane 
sections  remain  plane,  gives 

M  =  pl*t  .  I  =  h3/12  (5.2) 

in  which  I  is  the  second  moment  of  area,  h  the  wall  thickness,  and 
H  is  the  curvature.  If  the  cylinder  radius  is  a  and  the  radial  inward 
displacement  is  w(t>,  t) 


FIG.  5.2  IDEALIZED  STRESS,  STRAIN- 
RATE  LAW 


a  =  (w§9  wj/a^  (5.3) 

in  which  the  subscripts  9  denote  partial  differentiation  with  respect 
to  the  angular  coordinate  9.  The  final  inward  displacement  is  con¬ 
sidered  small  enough  to  allow  he  u  <e  of  the  initial  values  of  the  radius 
and  thickness. 

The  lack  of  uniformity  around  the  ring  of  the  mean  circum¬ 
ferential  stress  is  neglected.  Also,  it  is  assumed  that  the  slope  p 


of  the  stress  vs.  strain-rate  line  is  small  so  that  the  mean  circum¬ 
ferential  stress  may  be  considered  independent  of  time.  Consequently, 
a  constant  circumferential  force  5  is  taken  throughout  the  inward 
motion. 

5.  3  EQUATIONS  OF  MOTION 

From  Fig.  5.  3,  which  shows  an  element  of  the  ring  with  its 

attending  forces  and  moments,  the 
equations  of  equilibrium  or  motion 
(neg)ecting  rotary  inertia)  may  be 
written  as 

Q  =  (5.4) 

and 

Q\  +  Sc(\  =  “mwtt  (5- 

in  which  Q  is  the  shear  force,  dX 
the  arc  element  corresponding  to 
d6  but  now  subtending  an  angle  dcp, 
m  the  mass  per  unit  run  of  circumference,  and  the  subscripts  denote 
partial  differentiation.  is  the  curvature  and  is  given  by 

cp^  =  1/a  +  h  (5.  6) 

Combining  Eqs.  (5.  4)  and  (5.  6)  and  performing  the  operations 
necessary  to  eliminate  Q,  M  and  h  yields  the  following  equation  for 
w  ; 

^(ve000  +"99)t/«4  +  S[l/a  +  (w00+w)/a2]  =  -mwu  (5.7) 

It  is  convenient  to  introduce  the  dimensionless  displacement 
and  time  variables 

u  =  w/a  ,  t  =  plt/ma4  (5.8) 


FIG.  5.3  NOTATION.  ELEMENT 
OF  SHELL  IN  MOTION 
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It  is  an  approximation  to  the  duration  when  slight  wrinkling  occurs 
since  the  plastic  work  done  in  shortening  the  circumference  is  not 
influenced  to  the  first  order  by  the  slight  wrinkling. 

A  slightly  nonuniform  initial  velocity  may  be  represented  in 
the  form 

00 

w  (0,  o)  =  V  [  1  +  £  (a  cosn6  +  8  sinn0  )] 
con  n 

n-2 

»  (5.15) 

u  (0,  o)  =  V  [  1  +  T,  (a  cos  n0  +  8  sinn0)j 
T  o  n  n 

n=2 

Terms  with  n=l  are  omitted  since  they  do  not  contribute  to  the  defor¬ 
mation  of  the  circular  shape. 

The  solution  of  the  reduced  form  of  (5.  10)  may  be  taken  in  the 

form 

00 

u=  £  [fn(i)  cos  n0  +  gn(T)sinn0]  (5.16) 

n=2 

where  fR  must  satisfy  the  equation 

fn  +  n  (n  -  l)fn  -  (n  -  l)s  fn  =  0 

The  function  gn  must  satisfy  the  same  differential  equation. 

The  appropriate  solutions  for  f  and  g^  are 

\+  T  k" T 

f  =  B  e  n  +  C  e  n  (5.  17) 

n  n  n  '  ' 


Xn  T  O 

*n  =  V  +V  18> 
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where 


^  =  [n2(n2-l)/2][(l+4s^/^n‘*(n‘-l)h1/C-l] 


.4,2 


,1/2 


X‘  -  [n2(n2  -  l)/2]  [  1  +  4a2 
n 


!  |n4(n2  -  l)j 


(5.19) 


>1/2+l] 


The  constants  in  (5.  17)  and  (5,  18)  are  obtained  from  the  initial 
conditions  u(0,  o)  =  0  and  uT(0,  o)  given  by  (5.  15).  When  the  resulting 
functions  f  and  are  substituted  into  (5.  16)  and  the  resulting 
complementary  solution  added  to  the  particular  solution  (5.  13)  the  final 
solution  is 


u  =  -  1  +  cos  8 t  +  (v  / s  )sins  t  +  v  V  A  (t)(cx  cos  n0  +  B  sinnG)  (5.20) 
o  o  n  n  n 

n=2 


where 


\~  t  A  t 
A  (T)  =  (e  "  -« 


n 


(5.21) 


and  A*  being  given  by  (5.  19). 

From  a  comparison  of  the  expressions  (5.  15)  and  (5.20)  it  can 
be  seen  that  Aq(t)  may  be  considered  as  an  amplification  factor  which 
depends  on  the  harmonic  and  time.  Knowing  Aft(T)  determines  the 
effect  of  the  initial  velocity  perturbations  on  the  motion.  One  important 
property  of  the  amplification  factor  is  that  it  tends  to  zero  as  n  increases. 
In  fact  for  n  large  enough  Ar  ~  1/n*  .  This  ensures  the  convergence 
of  the  aeries  representation  (5.  20)  and  shows  that  the  higher  harmonics 
in  the  velocity  perturbation  will  have  little  influence  on  the  final  wrinkled 
state. 


5.5  STRAIN  REVERSAL 


The  foregoing  theory  is  valid  until  strain  reversal  begins. 
For  points  on  the  outside  and  inside  surfaces  the  circumferential 
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compressive  strain  is 

(1 -h/2a)(w/a -hWgg/2a^)  and  {1 +h/2a)(w/a +  hWgg/2a^) 

and  strain  reversal  begins  when  one  of  these  reaches  a  maximum 
positive  value  and  then  decreases.  When  this  occurs  the  strain  rate 
at  the  inner  or  outer  fiber  is  zero.  The  second  factors  in  the  above 
expressions,  in  dimensionless  terms,  are  u  +  hugg/2a  ,  and  from 
(5.20)  are  representable  in  the  form 


u  +  hUgg/2a  =  -  1  +  cos  st  (vq/s)  sin  st 


h  2. 


(5.22) 


+  v0  r  t  2aT  n  )  An^  *  ^anC°8  n9  +  9n  sinn6  ) 

n=2 


the  strain  history  is  thus  determined  by  (5.22). 


5.  6  RANDOM  VELOCITY  PERTURBATIONS 


From  (5.  15)  the  perturbation  of  the  initial  velocity  is 


u  (6,  o)  =  £  fa  cos  n0  +  b  sinnQ)  =  £  c  cos  (n8  -  0  )  (5. 23) 

t  n  n  n  n 

n=2  n=2 


where 


2  2  2  2 

a  =  v  a  ,  b  =  v  3  ,  and  c  =  v  (a  +  g  ) 

n  on  n  on  n  o'  n  n 


It  can  be  seen  from  (5.  20)  that  if  u^Q,  o)  is  given,  the  final  wrinkled 
shape  is  known.  With  a  statistical  approach,  Lindberg  [5.7,  5.8] 
demonstrated  that,  in  the  case  of  the  impact  buckling  of  thin  bars  and 
cylinders,  the  assumption  of  imperfections  in  the  initial  displacements 
of  the  white  noise  type  (displacements  composed  of  harmonics  of  equal 
amplitude,  random  phase  and  with  a  uniform  spectral  density)  led  to 
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a  good  description  of  experimentally  observed  wavelengths.  Here, 
a  similar  assumption  will  be  made  concerning  the  perturbation  of  the 
initial  velocity.  The  random  process  will  be  considered  stationary 
and  ergodic.  In  order  to  bring  out  the  analogy  between  the  present 
problem  and  the  relevant  parts  of  the  article  of  Rice  [  5.  9]  on  random 
noise  currents  the  following  change  of  variable  and  substitutions  are 
made  in  (5.  23' 

8  =  2rrcp/cp 


where 


0  £  cp  £  cp 


f  =  n/cp  and  <ju  =  2rrf 
n  n  n 


then  (5.  23)  becomes 


ut(8,  o)  =  £  cncos(uyp-  6n) 
n=2 


(5.24) 


Now  the  kinetic  energy  associated  with  the  perturbed  initial 
velocity  is 


2tt  $  00 

j-ma  j  u*(0)d6  =  y^nam)!  /  **  (cp)dcp  =  y(2nam)|-  £  c^ 
o  o  1 

2 

and  that  associated  with  the  n-th  component  is  proportional  to  cn/2  . 

Defining  w(f  )  as  a  kinetic  energy  density  (power  spectral 
n  2 

density)  it  follows  that  w(f  )Af  =  c  /2  where  &{  is  the  "bandwidth" 

7  n  n 

associated  with  the  n~th  component  of  u^,  so  that 

*f  =  fn+l  -  fn  =  1/(*> 
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For  a  white  noise  velocity  perturbation  c  =  c  ,  a  constant, 

n  2 

and  the  kinetic  energy  density  is  uniform  at  w(fR)  =  c  cp/2  .  Also  the 
phase  angles  6r  in  (5.24)  are  randomly  distributed  between  0  and 
2tt  . 


It  follows  from  (5.20)  that  the  perturbed  deformation  is 

N 

u(8,  t)  =  c  £  Ar(t)  cos  (u)n cp  -  6n)  (5.25) 

2 

4 

where  N  is  large  enough  to  consider  An  negligible  (An~l/n  for 
large  n  and  shows  a  behavior  analogous  to  a  filter). 


Applying  now  the  theorem  of  Rice  [5.9]  for  the  expected 
number  of  zeros  per  unit  angle  cp  and  converting  this  to  the  number 
of  wave 8  for  0  in  (0,  2rr )  gives 


— 1  1/2 


fn‘ 


A  d 
n  n 


UJ 

I 


AZd 
n  n 


(5.26) 


5.  7  NUMERICAL  ANALYSIS  AND  EXPERIMENTAL  RESULTS 

Figure  5.  4  is  a  stress -strain  diagram  obtained  from  a  static 
tensile  test  on  one  of  several  specimens  of  fully  annealed  1015  steel 
which  were  cut  from  cylinders.  It  demonstrates  the  lack  of  strain 
hardening  below  about  3  percent  strain.  According  to  the  experimental 
results  of  Manjoine  [5.  l]  the  material  is  extremely  sensitive  to  the 
rate  of  strain;  both  the  length  of  the  horizontal  or  ideally  plastic  portion 
of  the  stress -strain  curve  and  the  yield  stress  increase  with  strain 
rate.  Based  on  these  results  an  empirical  relation  between  yield  stress 
o  and  strain-rate  e(  has  been  suggested  by  Symonds  and  Bodner  (  5.  9] 
in  the  form 

o/oy  *  1  +  (et/D)1/p  (5.27) 
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FIG.  5.4  STATIC  STRESS-STRAIN  DIAGRAM 
FOR  FULLY  ANNEALED 
1015  STEEL 


and  the  value  of  0Q 
with  the  axis  e  =  0  . 


where  c  is  the  static  yield 

y 

stress,  and  p  and  D  are 
empirical  constants.  For  mild 

steel  these  are  taken  as  o  = 

2  7 
30,  000  lb /in  ,  p  =  5  ,  and 

D  =  40.  4  sec  *  . 

The  curve  (5.27)  is  shown 
in  Fig.  5.  5.  In  the  development 
of  the  theory  the  linear  strain - 
rate  law  (5.  1)  was  used  and  this 
is  here  taken  as  the  equation  of 
the  tangent  to  the  curve  (5.27) 
at  the  point  corresponding  to  the 
value  of  the  initial  strain-rate 
et(0).  An  example  of  such  a 
tangent  is  shown  in  Fig.  5.  5. 

The  slope  of  the  tangent  is 

H  =  oy/(PD1/pct1"1/p)  (5.28) 


is  given  by  the  point  of  intersection  of  the  tangent 


All  numerical  cases  correspond  to  experiments  that  arc 
described  later  and  the  main  results  are  listed  in  Table  5.  1.  The 
initial  velocities  VQ  are  obtained  from  impulse  calibration  experi¬ 
ments  lor  the  explosive  used  in  similar  configurations.  The  sequence 
of  operations  for  finding  the  predicted  number  of  waves  around  the 
cylinder  is  as  follows.  From  the  value  of  the  initial  strain-rate 
et(0)  =  V/a  is  found  and  hence  from  (5.27),  (5.28),  and  <5.  1).  and 
p  are  found.  Next,  the  dimensionless  constants  S  and  V  are 
determined  by  (5.  9)  and  (5.12).  Equation  (5.  14)  can  now  be  solved  for 
t  which  represents  the  duration  of  motion;  then  (5.8)  gives  the  actual 


approximate  duration.  The  spectrum  of  the  amplification  factor 
An(i)  follows  from  (5.21).  Finally,  the  expected  number  of  waves 
ir  computed  from  {5.  26). 


FIG.  5.5  STRAIN-RATE  LAWS 


Table  5.  1  contains  the  values  of  n  where  A  is  a  maximum 

n 

at  each  quarter  of  the  duration.  It  also  contains  the  values  of  at 
these  times.  This  is  done  to  show  the  variation  of  the  most  amplified 
harmonic  and  the  expected  number  of  waves  with  time.  In  many  c*3es 
the  variation  is  not  great.  Figures  5.  6  and  5.7  show  the  amplification 
spectra  for  the  cylinders  numbered  3a  and  4a.  For  case  4a,  it  can  be 
seen  from  Fig.  5.  7  th  -  there  is  a  rapid  growth  of  the  amplification 
factors  for  harmonics  about  the  value  n  =  12  and  that  above  n  =  30 
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virtually  no  amplification  occurs  (A^  —  1/n  for  large  n  ).  During 
the  final  three-quarters  of  the  duration  of  motion  the  values  of  n  at 
peak  move  from  n  =  8  to  n  =  12  while  the  expected  number  of 
waves  moves  only  from  =  12  to  N.  *  14  ,  The  magnification 
factor  curves  of  Figs.  5.  6  and  5.  7  are  quite  typical  and  similar  shifts 
in  the  most  amplified  harmonic  and  expected  number  of  waves  occur 
as  deformation  proceeds. 

Again  for  cases  3a  and  4a  random  initial  velocity  perturbations 
were  considered.  In  case  3a,  the  first  50  harmonics  of  constant  ampli¬ 
tude  c  and  random  phase  (white  noise)  were  employed  in  (5.23)  and 
the  representation  is  shown  in  Fig.  5.8(a).  The  resulting  perturbed 


FIG.  5.8  RESPONSE  OF  CYLINDER  3o  TO  WHITE  NOISE  VELOCITY  PERTURBATIONS 

(o)  Initial  velocity  perturbation,  (b)  Buckled  form  oftar  28  ;  tec,  (c)  Buci  1*4  form  oft#* 
42  fitec,  ofxl  (d)  Buckled  form  oftar  56  ptac 
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inward  displacements  in  the  form  u/c  were  obtained  from  (5,25)  at 
one -half,  three-quarters,  and  at  the  whole  of  the  duration  of  motion. 

They  are  shown  in  Figs.  5.8(b),  (c),  and  (d),  respectively  and  provide 
an  illustration  of  the  development  of  the  buckled  form.  The  higher 
harmonics  in  the  initial  velocity  perturbation  are  filtered  out  because 
of  the  form  of  the  spectrum  of  the  amplification  factor  (Fig.  5.  6).  In 
the  example  of  Fig.  5.8,  a  preferred  wavelength  exists  corresponding 
to  18  waves  around  the  cylinder,  the  expected  number  of  waves  if 

=  14,  and  one  experiment  gave  13  waves.  For  case  4a,  Fig.  5.7(a) 
shows  a  perturbed  initial  velocity  using  the  first  30  harmonics  of  equal 
amplitude  and  random  phase.  The  development  of  the  buckled  form  can 
be  seen  in  Figs.  5.7(b),  (c),  and  (d).  The  higher  harmonics  are 
filtered  out,  15  waves  develop  around  the  cylinder,  the  expected  number 
of  waves  is  =  14,  and  the  experimental  number  is  9  waves. 

The  amplitudes  of  the  initial  imperfections  are  considered  small 
enough  to  avoid  strain  reversal.  Experimental  values  for  the  number 
of  waves  comprise  the  right-hand  column  of  Table  5.  1  and  can  be  com¬ 
pared  with  the  expected  number  of  waves  of  the  most  amplified  harmonic, 
the  former  being  more  meaningful  but  the  latter  serving  as  an  indicator. 
The  degree  of  agreement  lends  credibility  to  the  postulated  mechanism, 
at  least  as  a  first-order  description. 

For  cylinders  3a  and  4a,  Figs.  5.  10  and  5.  11  show  the 
developed  final  shapes  and  can  be  compared  with  the  predicted  shapes 
of  Figs.  5.  8(d)  and  5.  9(d)  respectively.  Apart  from  the  lower  harmonics 
present  in  the  experiments  due  to  a  slight  lack  of  circularity  of  the 
cylinders  the  buckled  forms  are  similar.  For  these  two  cylinder  har¬ 
monic  analyses  (trigonometric  interpolation)  [5.  10]  of  the  experimental 
buckled  forms  were  carried  out  and  the  results  are  shown  in  Figs. 

5.  12  and  5.  13.  The  representation  was  taken  ii  the  form 

CO 

u(0,  t)  =  £  cncos(n6  -  cpn) 
o 
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with  data  recorded  at  every  degree  and  the  coefficients  c^  plotted 
against  n  .  The  large  coefficients  for  the  lower  harmonics  are  due  to 
the  initial  lack  of  circularity  of  the  cylinders.  In  case  3a  peaks  occur 
at  n  =  6,  9,  11,  and  15,  and  components  above  n  =  22  are  negligible. 
In  case  4a  peaks  occur  at  n  =  7,  10,  and  13,  and  components  above 
n  a  16  are  negligible.  No  major  conflict  arises  with  the  amplitude 
curves  of  Figs.  5.  6  and  5.  7. 


FIG.  5.9  RESPONSE  OF  CYLINDER  4a  TO  WHITE  NOISE  VELOCITY  PERTURBATIONS 
(a)  Initial  velocity  perturbation,  (b)  Buckled  form  after  28.5  /isec,  (c)  Buckled  form 
after  42.75  fxsec,  and  (d)  Buckled  form  after  57  jisec 
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5.  8  EXPERIMENTAL  TECHNIQUES 

The  buckled  cylinder  of  Fig.  5.  1  was  produced  with  the  experi¬ 
mental  arrangement  shown  in  Fig.  5.  14.  Enclosing  the  cylinder  is  an 
attenuator  of  1/ 4-inch -thick  neoprene  foam  and  a  layer  of  explosive. 

The  attenuator,  which  is  necessary  to  prevent  the  cylinder  from 
spalling,  extends  beyond  the  end  of  the  cylinder  to  provide  support 
for  the  explosive.  The  detonator  starts  from  an  electrical  signal  and 
is  connected  to  the  main  charge  by  several  strips  of  explosive.  Deto¬ 
nation  fronts  initiated  in  the  main  charge  by  the  strips  of  explosive 
expand  from  the  initiation  points  and  coalesce  into  a  single  ring-shaped 
detonation  front  in  the  "run-up"  part  of  the  main  charge.  The  detonation 
front,  which  is  of  the  order  of  a  few  mils  wide,  subsequently  sweeps 
over  the  cylinder  at  a  rate  of  23,  300  ft/sec. 


FIG.  5.11  BUCKLED  FORM  OF  CYLINDER  4o 
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DETONATOR 


FIG.  5.14  EXPERIMENTAL  ARRANGEMENT 


A  shock  wav*  in  induced  in  the  attenuator  by  the  detonation 
front  and  is  transmitted  from  the  attenuator  to  the  cylinder,  imparting 
an  impulsive  velocity  VQ  (ace  Table  5.  1). 

Maximum  variations  in  the  thickness  of  the  explosive  are  from 
21  to  5  percent  for  charges  from  12  to  5?  mils  thick,  the  larger  variations 
going  with  the  thinner  charges.  The  variations  in  thickness  are  gradual 
except  across  joints.  Variations  in  the  thickness  of  the  explosive  and 
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the  uncertainty  in  the  impulse  developed  by  the  explosive  combine  to 
give  an  uncertainty  of  the  order  of  10  percent  in  the  velocities  imparted 
to  the  cylinders. 

Table  5.  1  gives  the  properties  of  the  cylinders  tested  and  the 

initial  velocities  V  were  obtained  by  taking  an  impulse  constant  for 

°  5  3 

the  explosive -attenuator-target  configuration  of  2.  4  x  10  dyne-sec/cm  . 

5. 9  CONCLUSIONS 

In  its  present  form  the  theory  does  a  satisfactory  job  of  pre¬ 
dicting  the  buckled  forms  of  cylinders  of  rate -sensitive  material. 

There  appear  to  be  no  major  contradictions  to  the  postulated  mechanism 
of  deformation  although  it  would  be  desirable  to  have  more  experimental 
evidence  and  to  look  more  thoroughly  into  other  aspects  of  the  problem 
such  as  the  sensitivity  of  the  solution  to  the  strain-rate  law  and  to  the 
initial  velocity,  and  strain  reversal.  The  theory  as  it  stands  is  some¬ 
what  idealized  but  seems  adequately  compensated  by  simplicity  and 
usefulness. 
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6.  CORRELATION  OF  STRUCTURAL  RESPONSE  INFORMATION 


6.  1  INTRODUCTION 

Over  the  past  several  years  there  have  been  many  investigations 
of  the  failure  of  space -vehicle -type  structures  under  suddenly  applied 
external  surface  loads.  It  has  been  suggested  that  the  results  of  these 
investigations  be  correlated  to  maximize  their  usefulness  and  indicate 
areas  where  further  investigations  are  desirable. 

The  problem  of  correlating  the  results  of  structural  response 
investigations  has  at  least  two  aspects:  (1)  the  development  of  a  scheme 
for  comparing  the  results,  and  (2)  the  incorporation  of  the  results  into 
the  scheme.  The  following  presentation  is  based  on  this  division  of  the 
problem. 

6.  2  SCHEME  FOR  CORRELATING  THE  RESULTS  OF 

STRUCTURAL  RESPONSE  INVESTIGATIONS 

To  facilitate  the  correlation  of  the  results  of  structural  response 
investigations  it  is  highly  desirable  to  have  a  simple  and  direct  me.  rs 
of  comparison.  For  static  loads,  results  are  usually  summarised  by 
giving  the  load  magnitude  which  produces  a  critical  stress,  strain,  or 
displacement  in  the  structure.  The  same  format  appears  suitable  for 
dynamic  loads. 

Dynamic  loads  of  the  same  spacial  distribution  and  pulse  shape 
(time  variation)  which  produce  the  same  critical  stresses,  strains, 
and/or  displacements  in  a  structure  exhibit  a  simple  relationship  when 
plotted  in  the  amplitude -impulse  plane.  Figure  6.  1  illustrates  the 
relationship  between  amplitude  and  impulse  for  rectangular  pulses  which 
produce  the  same  maximum  displacement  of  a  linear  oscillator.  In  the 
region  where  the  curve  approaches  the  vertical  asymptote  the  load 
amplitude  becomes  unimportant  and  only  the  impulse  is  significant. 
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In  the  region  where  the  curve  approaches  the  horizontal  asymptote, 

the  impulse  becomes  unimportant  and  only  the  amplitude  is  significant. 

* 

In  the  intermediate  region  both  amplitude  and  impulse  are  significant. 

Figure  6.2  shows  the  amplitude -impulse  curves  for  rectangular 
and  triangular  pulses  which  produce  the  same  displacement  of  a  single¬ 
degree -of -freedom  system  of  a  rigid-plastic  material.  The  curves 
have  the  same  general  shape  as  that  of  Fig.  6.  1. 

Tests  of  more  complex  systems  such  as  cylindrical  shells, 
which  have  a  variation  in  response  mode  with  load  amplitude  and 
impulse  (for  loads  of  the  same  spacial  distribution  and  pulse  shape), 
exhibit  the  same  general  relationship  as  that  of  Figs.  6.  1  and  6.2 
between  the  amplituae  and  impulse  combinations  which  produce  the 
same  parmanent  displacement. 

The  response  modes  of  the  preceding  examples  are  not  sensitive 
to  stress  waves  and  are  termed  unitary  modes.  Leads  which  excite 
modes  which  are  sensitive  to  stress  waves,  termed  laminar  modes, 
are  also  conveniently  displayed  in  the  amplitude -impulse  plane. 


For  a  linear  oscillator,  the  information  given  by  the  amplitude-impulse 
curve  (Fig.  6.  1)  is  related  to  that  given  by  the  shock  spectrum  (6.  1  ] . 
For  a  load  of  given  shape,  the  shock  spectrum  gives  the  response  of 
a  linear  oscillator  as  a  function  of  the  frequency  of  the  oscillator  for 
a  load  of  unit  amplitude,  while  the  amplitude -impulse  diagram  gives 
the  amplitude  and  impulse  combinations  which  evoke  a  given  response 
from  an  oscillator  of  given  frequency. 

In  a  complex  structure  many  modes  are  excited  and  a  one-to-one 
correspondence  between  the  amplitude-impulse  diagram  and  the  shock 
;p~ctrum  does  not  exist.  For  such  structures,  the  amplitude -impulse 
diagram  displays  the  loads  which  evoke  a  ^ivtn  peak  response  in  the 
structure  ,  including  effects  of  all  modes  of  deformation,  even  if  the 
response  is  nonlinear  and  linear  superpo.c  ition  is  invalid.  The  shock 
spectrum  gives  the  peak  response  of  each  mode,  and  in  some  cases 
(the  response  must  be  linear  and  the  frequency  spread  must  be  limited) 
these  can  be  superposed  to  determine  an  approximate  peak  response 
,  of  the  structure. 
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FIG.  6  2  AMPLITUDE-IMPULSE  DIAGRAM  FOR  A  ON  E-DEGREE - 

OF -FREEDOM  SYSTEM  OF  A  RIGID  PERFECTLY -PLASTIC 
MATERIAL.  Hw*  o  ia  yialS  strait,  K  it  wt>H  thickneti,  o  it 
cyMnd«r  •vdlut,  e  it  itroir,  and  p  it  dtntiiy. 


Figure  6.  3  shows  schematically  the  amplitude -impulse  diagram  for  the 
loads  required  to  produce  a  given  deformation  o'  a  thin  cylindrical  shell 
contained  *ithin  a  close-fitting  sleeve  which  ia  loaded  uniformly  on  the 
exterior.  The  impulse  is  constant  (vertical  part  of  curve)  until  the  load 
duration  exceeds  twice  the  transit  time  through  the  inner  shell.  For 
monotonicaily  decreasing  pulses,  »t  longer  durations  the  amplitude 
decreases  and  approaches  a  limiting  value. 

Figure  6.4  shows  the  amplitude -impulse  curve  for  uniform 
loads  of  rectangular  shape  required  to  produce  cracking  (from  hoop 
tension)  of  a  thin  cylindrical  shell  by  rebound  from  an  internal  mandrel. 
The  jagged  nature  of  the  cunt  is  due  to  peculiarities  associated  with 
the  rapid  decay  of  rectangular  pulses. 
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fiC-  o.3  SCHEMATIC  AMPLITUDE-IMPULSE  DIAGRAM  FOR  A  THIN 
CYLINDRICAL  SHELL  CONTAINED  WITHIN  A  CLOSE- 
FITTING  SLEEVE  LOADED  ON  THE  FXTERIGR 

fVom  the  foregoing,  it  appears  that  amplitude -impulse  diagrams 
afford  a  convenient  basis  lor  comparison  of  the  results  of  structural 
response  investigations  having  the  ultimate  objective  of  predicting 
failure  loads.  A  particularly  desirable  feature  is  that  the  asymptotes 
of  the  amplitude -impulse  curves,  which  correspond  to  short -duration 
and  long-dura'ion  loads,  are  generally  the  easiest  part  of  the  curves 
to  establish  theoretically  and  experimentally.  Also,  for  most  structures, 
the  general  shape  of  the  curve  will  be  Similar  to  that  of  Figs.  6.  1  to 
n.4.  and  hence  interpolation  is  simple.  Another  desirable  feature  is 
that  amplitude  and  impulse  are  easily  related  to  weapon  yield  and 
range. 


With  this  approach,  to  completely  describe  the  failure  loads 
for  a  given  structure  with  a  given  failure  criteria  require*  a  set  of 
amplitude -impulse  curves,  one  curve  for  each  load  distribution  and 
pulse  shape.  Such  a  set  of  curves  would  appear  to  satisfy  the  needs  of 
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FIG.  6.4  SCHEMATIC  AMPLITUDE-IMPULSE  DIAGRAM  FOR  RECTANGULAR  PULSES 

TO  PRODUCE  CRACKING  (from  hoop  tension)  OF  A  THIN  CYLINDRICAL  SHELL 
BY  REBOUND  FROM  AN  INTERNAL  MANDREL 


the  designer,  who  must  design  structures  to  carry  certain  loads,  and 
of  the  analyst,  who  must  determine  the  failure  loads  of  given  struc¬ 
tures  . 


6.  3  CORRELATION  OF  THF,  RESULTS  OF  STRUCTURAL 
RESPONSE  INVESTIGATIONS 


Of  t!  vast  amount  of  structural  response  information  avail¬ 
able,  only  a  very  small  part  w  relevant  to  the  problem  of  determining 
failure  loads  ox  s pace -vehicle -type  structures  under  suddenly  applied 
external  surface  loads.  Most  of  the  available  structural  response 
information  is  for  linear  elastic  systems  and  is  of  secondary  importance 
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for  the  present  application  since,  for  the  structures  of  greatest  interest, 
failure  generally  involves  significant  permanent  deformation.  More¬ 
over,  of  the  available  information  on  permanent  deformation  of  struc¬ 
tures,  most  is  for  static  loads. 

Information  on  permanent  deformation  of  space -vehicle -type 
structures  under  suddenly  applied  external  surface  loads  appears  to 
stem  from  three  main  rces:  (1)  investigations  of  effects  of  air 
blasts  on  airplane  structures,  {2}  investigations  of  the  effects  of  under¬ 
water  blasts  on  ship  and  submarine  hulls,  and  (3)  recent  investigations 
on  the  effects  of  suddenly  applied  external  surface  leads  on  space - 
vehicle -type  structures  The  re'  ults  of  the  investigations  on  aircraft 
are  of  some  interest  for  lightweight  sat-.  Mtes,  but  are  of  very  limited 
use  for  the  relatively  thick-walled  structures  of  re-entry  vehicles. 

The  results  of  investigations  on  ship  and  submarine  hulls  should  provide 
significant  information  on  re-entry  vehicle  structures.  Unfortunately, 
this  information  has  been  inaccessible.  Therefore,  the  following 
discussion  is  based  on  the  recent  investigations  undertaken  specifically 
to  obtain  additional  information  on  the  response  of  space -vehicle -type 
structures  under  suddenly  applied  external  surface  loads. 

The  basic  structural  element  of  space  vehicles  in  general  is 
the  cylindrical  shell,  or  conical  shell  of  small  angle.  The  first  known 
attempt  to  determine  bounds  on  a  wide  range  of  suddenly  applied 
external  surface  loads  which  produce  a  given  permanent  deformation 
of  a  cylindrical  shell  is  that  described  in  [6.2],  and  was  concerned 
with  smoothly  varying  load  distributions  The  bounding  curves  found 
are  shown  in  Fig.  6.  5.  The  upper  bound,  shown  for  rectangular  and 
triangular  pulse  shapes,  are  the  loads  which  produce  10  percent  defor¬ 
mation  of  a  uniformly  loaded  cylindrical  shell  of  a  rigid-plastic  material 
which  responds  in  uniform  radial  motion  (i.e.,  buckling  is  excluded). 
Since  in  this  mode  the  maximum  amount  of  energy  is  absorbed  in  plastic 
work  per  unit  deformation,  and  no  actual  load  could  produce  such  an 
idealized  response,  this  is  an  absolute  upper  bound  on  actual  loads 
which  produce  the  same  deformation,  for  any  cylindrical  shell. 
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FIG.  6.5  BOUNDING  CURVES  ON  SUDDENLY-APPLIED  EXTERNAL 
SURFACE  LOADS  OF  SMOOTHLY  VARYING  DISTRIBUTION 
TO  PRODUCE  10  PERCENT  DEFORMATION  OF  A 
CYLINDRICAL  SHELL.  ($••  Caption  to  Fig.  6.1  for  definition*.) 

The  lower  bound  la  baaed  on  limited  experimental  results 
[  6.  3] ,  [6.6],  and  in  view  of  the  extensive  effort  being  made  to  obtain 
additional  information  on  failure  loads  of  cylindrical  shells,  should  be 
regarded  as  tentative. 

The  effort  alluded  to  above  is  that  being  undertaken  at  SRI 
under  Contract  P.  O.  24-14517  under  AF  \694)-655  with  Lockheed 
ivlissiles  and  Space  Company  for  the  Air  Force.  This  is  a  combined 
experimental  and  theoretical  effort  with  the  objective  of  establishing 
the  amplitude -impulse  curves  for  suddenly  applied  external  surface 
loads  of  smoothly  varying  distribution  which  produce  10  percent  defor¬ 
mation  of  monocoque  cylindrical  shells  and  multilayer  cylindrical 
shells  typical  of  space  vehicles.  The  experiments  thus  far  have  shown 
that  buckling  is  a  predominant  failure  mechanism  over  the  entire  load 
range.  The  theory  which  is  being  developed  is  based  on  this  mechanism. 
A  report  on  this  work  is  expected  to  be  available  in  late  19Cr 
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Except  for  the  extensive  experimental  effort  by  BRL  [  6.  6] ,  the 
results  of  which  were  included  in  determining  the  lower  bound  shown  in 
Fig.  6.5,  the  results  of  other  investigations  on  cylinders  have  not  yet 
been  cast  in  the  form  of  amplitude -impulse  diagrams. 

Other  structural  elements  of  interest  for  space  vehicles  are 
spherical  shells,  beams,  and  plates.  For  spherical  caps,  an  upper 
bound  can  be  obtained  on  the  same  basis  as  that  for  cylindrical  shells. 

The  resulting  curve  is  given  by  Fig.  6.5  with  pQ  =  2oh/a  and  h 

No  data  or  theory  are  known  which  would  establish  a  curve  of  load  ampli¬ 
tude  and  impulse  to  produce  a  given  permanent  deformation  of  spherical 
shells. 


In  a  recently  submitted  report  [6.4],  theoretical  ampli.ude- 
impulse  curves  are  presented  for  rigid-plastic  beams  and  plates  under 
a  -variety  of  boundary  conditions  and  loads.  However,  the  curves  depart 
by  factors  of  3  to  10  from  experimental  results.  Considerable  more 
work  is  needed  to  refine  the  amplitude -impulse  curves  for  beams  and 
plates . 


Amplitude -impulse  curves  for  failure  (in  hoop  tension)  of  a  thin 
cylindrical  shell  by  rebound  from  an  internal  sleeve  for  uniformly  distri¬ 
buted  pulses  of  rectangular  form  are  given  in  [  6.  5].  This  work  will  be 
extended  in  the  near  future  for  uniformly  distributed  pulses  of  triangular 
and  exponential  forms. 

6. 4  CONCLUSIONS 

Amplitude  impulse  diagrams  appear  to  offer  a  meaningful  and 
convenient  way  to  represent  dynamic  failure  loads.  If  this  scheme  is 

adopted,  amplitude -impulse  curves  should  be  established  for  failure  > 

i 

loads  for  all  important  structural  elements  of  space  vehicles. 

Presently  available  theoretical  and  experimental  information 
for  establishing  failure  loads  of  space-vehicle -type  structures  is 
severely  limited.  Investigations  now  underway  should  significantly 
improve  this  situation. 
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7.  OTHER  INVESTIGATIONS 


7.  1  INTERACTION  OF  SHOCK  EFFECTS  AND  STRUCTURAL 

EFFECTS 

When  a  shell  structure  is  exposed  to  an  external  impulsive 
surface  load  a  shock  wave  is  generated  which  propagates  in  from  the 
exposed  surface.  If  the  load  is  sufficiently  intense,  the  structure  will 
be  destroyed. 

Two  basic  types  of  damage  from  impulsive  surface  loads  have 
been  identified;  shock  damage  and  structural  damage.  Shock  damage 
is  that  due  to  stress  wave  interactions  and  generally  occurs  on  the  first 
or  second  transit  of  the  stress  wave  through  the  wall  of  the  structure, 
and,  for  common  space  vehicles,  this  is  usually  within  a  few  micro¬ 
seconds  after  the  application  of  the  load.  Structural  damage  is  that  due 
to  the  average  velocity  of  the  shell  of  the  structure  and  generally  occurs 
50  psec  or  more  after  the  application  of  the  load,  for  common  space 
vehicle  structures. 

In  investigating  the  response  of  space  vehicle  structures  to 
impulsive  surface  loads,  it  has  become  common  practice  to  study  shock 
efiecis  and  structural  effects  separately.  This  division  is  certainly 
valid  for  shock  effects,  which  occur  too  early  to  be  influenced  by  struc¬ 
tural  action.  However,  shock  effects  may  influence  structural  effects. 

The  division  of  the  response  problem  into  shock  effects  and 
structural  effects  is  based  on  the  vastiy  different  time  scales  of  the  two 
phenomena.  Except  for  obvious  interaction  effects  such  as  spalling, 
cracking,  etc.,  this  division  assumes  that  the  stress  waves  have  sub- 
tided  to  insignificant  amplitudes  before  structural  deformation  becomes 
significant. 


For  two-layered  structures  of  nonmetallic  heat  shield  and  metal 
base  structure,  it  has  been  shown  [7.  l]  that  for  elastic  stresses, 
about  five  wave  transits  through  the  layers  is  sufficient  to  attenuate  the 
stress  amplitude  of  the  stress  waves  to  about  20  percent  or  less  of  the 
initial  amplitude.  For  common  space  vehicles,  this  requires  less 
than  50  psec.  Hence  for  such  structures  the  separation  of  the  response 
into  shock  effects  and  structural  effects  is  valid. 

The  attenuation  of  stress  waves  in  single-layer  structures 
depends  on  details  of  the  equation  of  state  of  the  material.  This  was 
not  investigated  on  the  present  contract. 

7.  2  PREDICTION  OF  ACCELERATIONS  FROM  SHOCK  LOADING 

Accelerations  from  suddenly  applied  surface  loads  can  produce 
accelerations  which  damage  internal  components  of  space  vehicles.  It 
would  be  highly  desirable  to  have  a  means  of  predicting  the  accelerations 
from  a  given  load.  However,  experiments  on  a  mocked-up  satellite 
vehicle  undertaken  on  Contract  AF  29(fc01)«6435  indicate  that  internal 
accelerations  are  very  sensitive  to  the  local  construction,  so  much  so 
that  it  appears  unlikely  that  a  useful  general  method  of  predicting 
accelerations  even  for  simple  structures  can  be  expected. 
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